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ABSTRACT 

We present new continuum and line observations, along with modelling, of the faint (6-8)Myr old TTauri star ETCha belonging 
to the 7] Chamaeleontis cluster. We have acquired Herschel/Pacs photometric fluxes at 70yum and 160yum, as well as a detection of 
the [OI] 63yum fine-structure line in emission, and derived upper limits for some other far-IR OI, CII, CO and 0-H2O lines. These 
observations were carried out in the frame of the Open Time Key Programme GASPS, where ET Cha was selected as one of the science 
demonstration phase targets. The Herschel data is complemented by new simultaneous Andicam B-K photometry, new Hst/Cos and 
Hst/Stis UV-observations, a non-detection of CO J = 3 — > 2 with Apex, re-analysis of a Ucles high-resolution optical spectrum 
showing forbidden emission lines like [OI] 6300 A, [SII] 6731 A and 6716 A, and [Nil] 6583 A, and a compilation of existing broad- 
band photometric data. We used the thermo-chemical disk code ProDiMo and the Monte-Carlo radiative transfer code MCFOST 
to model the protoplanetary disk around ETCha. The paper also introduces a number of physical improvements to the ProDiMo 
disk modelling code concerning the treatment of PAH ionisation balance and heating, the heating by exothermic chemical reactions, 
and several non-thermal pumping mechanisms for selected gas emission lines. By applying an evolutionary strategy to minimise the 
deviations between model predictions and observations, we find a variety of united gas and dust models that simultaneously fit all 
observed line and continuum fluxes about equally well. Based on these models we can determine the disk dust mass with confidence, 
Afdust ~ (2 - 5) x 10~ 8 M whereas the total disk gas mass is found to be only little constrained, M gas a (5 x 10~ 5 - 3 x 10 -3 ) M Q . Both 
mass estimates are substantially lower than previously reported. In the models, the disk extends from 0.022 AU (just outside of the 
co-rotation radius) to only about 10 AU, remarkably small for single stars, whereas larger disks are found to be inconsistent with the 
CO J = 3 — » 2 non-detection. The low velocity component of the [OI] 6300 A emission line is centred on the stellar systematic velocity, 
and is consistent with being emitted from the inner disk. The model is also consistent with the line flux of H 2 v=l— > S(l) at 2.122yum 
and with the [OI] 63 /jm line as seen with Herschel/Pacs. An additional high-velocity component of the [OI] 6300 A emission line, 
however, points to the existence of an additional jet/outflow of low velocity 40-65 km/s with mass loss rate ss 10~ 9 M Q /yr. In relation 
to our low estimations of the disk mass, such a mass loss rate suggests a disk lifetime of only ~ 0.05 - 3 Myr, substantially shorter 
than the cluster age. If a generic gas/dust ratio of 100 was assumed, the disk lifetime would be even shorter, only ~ 3000 yrs. The 
evolutionary state of this unusual protoplanetary disk is discussed. 

Key words. Stars: pre-main sequence; Protoplanetary disks; Astrochemistry; Radiative transfer; Line: formation; stars: individual: 
ETCha 



1. Introduction 



Gas-rich dust disks around young stars (hereafter, protoplane- 
tary disks) provide the raw material to build up new planets. 
The physical, thermal, and chemical conditions in the disk, the 
timescale over which the gas disperses, and the physical mecha- 



nisms contributing to the gas dispersal are keys to understanding 
what type of planets can form and on what timescales. 

Significant progress has been made in the past few years in 
measuring the dispersal timescale of the dust component of the 
disk. Infrared surveys of nearby star-forming regions and associ- 
ations have established that the frequency of o ptically thick dus t 
disks decreases exponentially with time (e.g. iMamaiekL I2009). 
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By an age of lOMyr only a few percent of pre-main sequence 
Sun-like stars ( T Tauri stars) stil l retain an optically thick dust 
disk (see e.g. iHernandez et al.L l2008t iPascucci & Tachibanal 
2010, for reviews). Since the near-mid infrared excess (A 5> 
30 ^m) is sensitive to the presence of small dust grains, not 
larger than a few microns in size, these observations effectively 
trace the dispersal of small grains within about 10AU from T 
Tauri stars. Millimetre observations, tracing colder dust at hun- 
dreds of AU from the central star, indicate a similarly fast clear- 
ing f or the outer disk, wit hin about 10-30Myr for T Tauri 
stars dCarpenter et al.ll2005l) . There is growing observational ev- 
idence that the dust disk lifetime depends on stellar mass. Disks 
around intermediate-mass (k, 1.5 M ) stars disperse in less than 
lOMyr, whereas disks around low-mas s stars (M dwarfs and 
brown dwarfs) persist for longer times dCarpenter et al.l 2006; 
ICurrie et all l2007t iRiaz & GizisL l2008h . 

Due to observational challenges in detecting gas lines from 
disks and difficulties in interpreting them, much less is known 
about the evolution of the gas component of the disks. Three ob- 
servables point to a dispersal timescale similar to (or possibly 
shorter than) the dust dispersal timescale: the ex ponential de- 
creas e with time in the frequency of accreting stars dFedele et all 
2010); the non-detections of infrared gas lin es from abundant 
molecules and atoms in ten uous dust disks dHollenbach et all 
12005b IPascucci et all H006); upper limits on the H2/dust mass 
ratio of less than 10 in two ~ 12 Myr old edge-o n disks 
dLecavelier des Etangs et~alll2001blRoberge et all 120051 ). 

The aim of the Herschel Open Time Key Program "Gas in 
Protoplanetary Systems" (GASPS, Dent et al. in prep.) is to pro- 
vide new insights into the chemical and gas temperature struc- 
ture of protoplanetary disks, the gas/dust ratio, the gas disper- 
sal timescale, and disk evolution. GASPS will acquire a large 
sample of sensitive far-infrared Herschel/Pacs spectra for 240 
disks in nearby star-forming regions and associations that span 
the critical l-30Myr age range over which disks are known 
to disperse. The primary signatures of the gas in the disk are 
expected to be the forbidden [OI]63.2//m, [OI] 145.5 ^m, and 
[CII] 157.7yum lines, as well as some CO and H2O lines. The 
first GASPS papers have shown that (1) the [OI] 63 /mi line can 
be used as primary g as indicator and is of ten detected toward 
protoplanetary disks (Mathe ws et all 1201 Ol) . (2) a combination 
of far-IR and (sub-)millimetre gas lines provides a promis- 
ing tool to estima te the total gas mass of protoplanetary disks 
(Pinte et all l2010). and (3) detailed models of individual sources 
allow to characterise the disk structure and shape, and the dust 
and gas compone nts of protoplanetary disks (Me eus et al 
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lThietalll2010ah . 

In this paper, we present an analysis of the circumstellar disk 
of ET Cha, an approximately 8 Myr old late-type T Tauri star, 
with the goals of characterising in detail its dust and gas con- 
tent. ETCha is one of the few nearby re latively old stars still 
posse ssing an optically thick du sk disk dSicilia-Aguilar et all 
120091) and still accreting disk gas dLawson et alll2004l) . TW Hva 
and PDS 66 are two other well-known old stars with proper- 
ties similar to ETCha. Both disks have been studied in de- 
tail and show evidence of evolution with respect to 1 - 2 Myr 
old T Tauri disks in Taurus, for example, depleted inner disk 
in TW Hva dCalvetet all 120021) and flatter disk structure for 
PDS 66 dCortes et all 12009). Both disks are likely to have too 
low disk masses to form giant planets at this evolutionary stage. 
ETCha would be the third such old disk system where obser- 
vational data allows for an in-depth-study of its dust and disk 
properties. 



The paper is structured as follows. Section 2 provides 
an overview of the prior knowledge of the source. We then 
present new multi- wavelength observations of ET Cha in Sect. 3. 
Section 4 presents a detailed dust and gas disk model for ET Cha. 
We describe the main results of our models in Sect. 5. Finally, 
we discuss some critical aspects of the modelling, and the im- 
plications of both models, in Sect. 6, before we finish the paper 
with our conclusions in Sect. 7. 



2. ETCha: an old TTauri star with active accretion 

ETCha (2 MASS J0843 1857-7905 181, ECHA J0843. 3-7905, 
also referred to as RECX15B is a low-mass T Tauri star that was 
identified by dLawson et all 120021) as a membe r of the nearby, 
8 Myr old rj Chamaeleontis moving group (Mamaiek et al., 
2000fl The associa t ion is located only 97 pc away from the 
Sun (Mamaiek et al., 1999) and is virtually unaffected by ex- 
tinction (Luhman& Steeghs. 120041) . an ideal set of conditions 
to study circumstellar disks in detail. A slighly smal l er dis tance 
to ETCha of 94.3 pc was reported by Ivan Leeuwenl (12007). but 
we have used the earlier and better known value of 97 pc for the 
modelling in this paper. Brandek er et al.l (12006 ) obtained high- 
angular resolution images of ET Cha and concluded that it has 
no companions (brown dwarfs) outside of 10 AU (30 AU). 

ET Cha is one of the few association members that possess 
a circumstellar disk, as indicated by a series of Spitzer obser- 
vations that revealed the presence of dust including strong mid- 
infrared silicate features (Megeat h et all 1200 5: Bouw man et all 
12006b iGautier et all [20081: ISicilia-Aguilar et all 120091) . Despite 
the age of the system, the infrared colours of the source are remi- 
niscent of those of much younger (1-2 Myr) circumstellar disks. 
Furthermore, optical spectroscopy of ET Cha has shown that it is 
undo ubtedly accret i ng, wi t h a very st r ong and broad H g emission 
line dLawson et all 12002b iLvo et all 12004b iLuhman & Steeghs, 



2004). Based on the observed H ff line, the mass accretion rate 
of ET Cha has been estimated to be 10~ 9 M G /yr dLawson et all 
2004), and the disk inclination (by modelling the H a line pro- 
file) to be about 60° as measured from face-on. The spectra also 
reveal a series of forbidden optical emission lines ([OI], [SII], 
[Nil]) that unambiguously indicate the presence of a jet/outflow. 
In addition, the stellar absorption lines in these spectra allowed 
for precise spectral typing of the central star; all estimates agree 
with a M3-M3.5 spectral classification. 

Among all members of the 77 Cha association, ET Cha shows 
the largest variability in the visible of order 0.3-0.4 mag, which 
places it at the high end of the var iability distribution of WTTS 
(see Fig. 1 in lGrankin et al.. 2008). The main feature in the ob- 
served lightcurve is a ~ 12 day period which has been found in 
two consecutive years. lLawson et al.l (120021) also noted a flare 
lasting for about 1.7 days. The more regular 12 day variations 
are most easily interpreted in terms of an accretion hotspot co- 
rotating with the star. However, our re-analysis of optical absorp- 
tion line profiles (see Sect. 13. 6t suggests that the stellar rotational 
period is much shorter, around 2 days. Therefore, the physical 
origin of the ~ 12 day period remains uncertain. The stellar vari- 
ability casts some doubt on the derivati on of stellar param eters, 
because the photometry reported by lLawson et al.l (|2002) was 
taken "near maximum light", when e. g. an accretion hot spot 



1 The ROSAT survey reported by (Mamaiek ~et"aTl Il999t l200(ih lists 
only RECX1-12. Lables 13-15 have been used to denote three post- 
ROSAT stars discovered in or near the cluster core, including ET Cha. 

2 We note that Luhman & Steeghs ( 2004) derived an age of the 77 Cha 
association of only 6+' Myr. 
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Table 1. Observed Line Fluxes [10 18 W/m 2 ] with Herschel/Pacs and Apex. Detection are listed as Fi + cr whereas non-detections 
are listed as <3cr. 



[01] 


[01] 


[CII] 


o-H 2 


o-H 2 


o-H 2 


CO/ = 36^35 


CO/ = 33^32 


CO/=29^28 


CO/=3^2 


63.18 pm 


145.52pm 


157.74 pm 


78.74 pm 


179.53 pm 


1 80.49 pm 


72.84 pm 


79.36 pm 


90.16pm 


866.96 pm 


30.5 + 3.2 


<6.0 


<9.0 


<30 


<5.0 


<5.2 


<8.0 


<24 


<9.6 


<0.05 



r|Chal5 [01] 63nm 

1.4 I 1 1 1 1 1 1 




-0.6 1 1 1 1 1 1 1 1 1 1 1 — ! 1 

62.9 62.95 63 63.05 63.1 63.15 63.2 63.25 63.3 63.35 63.4 63.45 

X (|im) 

Fig. 1. Herschel/Pacs line observations, with overplotted 
Gaussian to measure the integrated line flux of [OI] 63.18pm. 
All other line observations are non-detections as listed 
in Table Q] Errorbars indicate the 1 cr noise level. The 
[OI] 63.1 8 pm line profile is dominated by the spectral resolu- 
tion of the PACS spectrometer R » 3000. 

could contribute significantly to the observed flux. Furthermore, 
it prevents using optical and near-infrared colours in case of non- 
simultaneous data. For this reason, we have obtained new, simul- 
tan eous photometry, which is present ed in Sect. 13. 21 

iRamsav Howat & Greaves! d2007l) detected a clear gas signa- 
ture from the disk of ETCha in form of the H2 v=l— > S(l) line 
in emission at 2. 122 yum with Gemini/Phoenix, with an integrated 
line flux of (2.5 ±0.1) x 10 -18 W/m 2 . These high angular resolu- 
tion observations showed a narrow line (FWHM = 1 8 + 1 .2 km/s) 
centred on the stellar velocity to within ~ 1 km/s. No angu- 
lar offset between the li ne and the star was detected a t the 
level of 4 AU. Therefore, Ramsav Howat & Greavesl (120071) ar- 
gue tha t this line is emitted by H? gas in Keplerian rotation at 
~2 AU. lRamsav Howat & Greavesl d2007l) observed 3 other disk- 
bearing members of the 77 Cha associati on, but ET Cha wa s the 
only one with detectable H? em ission. iBarv et alJ (120081) and 
iMartin-Zaidi et alJ {2009, 2010) reported on several detections 
of H2-lines toward other sources where the emission is also 
likely originating from the disk rather than from an outflow. 

3. Observations and Data Reduction 

3.1. Herschel/PACS 

Herschel/Pacs observations were obtained for ETCha during 
the Science Demonstration Phase. The photometric observations 
were obtained in the scan map mode in the blue (70 yum) and 
the red (160 pm) filters. Two different scan map angles were 
used, 45° (obsid 1342187338, 133s) and 63° (obsid 1342189366, 
220s). Both scans were obtained at a medium scan speed of 
20"/s, with a cross scan step of 5" and scan leg length of 



3'. The number of scan legs was 8 for the 63° scan, and 
4 for the 45° scan. For spectroscopic observations, a 1669s 
PacsLineSpec (obsid 1342186314) and a 5150s PacsRangeSpec 
(obsid 1342187019) were obtained. The PacsLineSpec pro- 
vides two simultaneous spectra at wavelengths 62.93-64.43 pm, 
and 180.76 - 190.29 pm. PacsRangeSpec covers six spectral 
ranges of 71.81 -73.28 pm, 78.37-79.76 pm, 89.28-90.48 pm, 
143.59- 146.53 pm, 156.70- 159.47pm, and 178.51 -180.96pm. 
Observations were taken in the chop-nod mode, with a small 
2" dither. The target was centred at the central spaxel of the 
9.4" x 9.4" grid of the Pacs Integral Field Unit. The data was 
reduced using the Herschel Interactive Processing Environment 
(Hipe; Ott 2010) developer build version 3.0.1212, and the data 
reduction scripts provided at the Herschel data reduction work- 
shop held in January 2010. 

For the photometric data, a mosaic was created from the two 
scan maps. Aperture photometry was performed using an aper- 
ture radius of 16"in the blue, and 19.2"in the red. An aperture 
correction of 0.922 (blue), and 0.885 (red) was applied to the 
photometry. The aperture corrections were obtained from the 
Pacs PhotChopNod Release Note (Feb. 22, 2010). The flux cal- 
ibration uncertainty is estimated to be 5% in the blue and 10% 
in the red. For the spectroscopic data, we extracted the spectrum 
from the central spaxel, and then applied an aperture correction 
in order to minimise the flux loss in the neighbouring spaxels. 
Spectra from the central spaxel were extracted for both the A 
and the B nods. We then applied wavelength binning to each nod 
spectrum, using a bin size that is half the width of the instrumen- 
tal resolution. The final spectrum is the mean of the wavelength- 
binned spectra from the two nods. The absolute flux calibra- 
tion uncertainty is estimated to be 40%. We have detected the 
[OI] 63.2 pm emission line for ET Cha, while all other lines are 
undetected (see TableQ}. We used the IDL routine MPFITPEAK 
to fit an error-weighted Gaussian to the observed [OI] 63.2 pm 
line, and measured the integrated flux of the Gaussian line fit. 
The 1 -cr error to the line flux was calculated by setting the height 
of the Gaussian equal to the continuum rms value, and the width 
equal to the instrumental resolution. The continuum emission at 
the rest wavelength of 63.18 pm was estimated by fitting a first- 
order polynomial to the spectral region. 

3.2. CTIO/ANDICAM photometry 

We obtained new simultaneous optical/NIR photometry using 
the dual-channel ANDICAM instrument on the CTIO 1.3 m 
telescope on March 6, 2009. Photometric calibration was per- 
formed using the PG 1047 Landolt field. Both the optical and 
near-infrared data were reduced using standard procedures (flat- 
fielding, cosmetic cleaning, shift-and add). The photometry on 
ETCha and the photometric standard was extracted within a 
3"aperture and airmass and colour corrections were applying 
using coe fficient s from the ANDICAM websit^l in the optical 
and from Frogel (1998) in the near-infrared. The resulting pho- 
tometry is listed in Table|2] along with the mid- and far-infrared 



3 http://www.astro.yale.edu/smarts/smarts 1 3m/photometry.html 
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Table 2. Photometric data of ETCha. 





mag. 


F v [Jy] 


instrument 


ref. 


used data ... 


0.091 -0.111 


- 


9.81e-6 


Fuse (scaled) 


GH 


0.111 -0.145 


- 


2.85e-4 


Hst/Cos 


GH 


0.145-0.205 


- 


1.69e-4 


Hst/Cos/Stis 


GH 


0.442 (B) 


15.64 


0.00195 ± 0.0001 


Ctio/Andicam 


GD 


0.55 (V) 


14.68 


0.0046 ± 0.00025 


Ctio/Andicam 


GD 


0.66 (R) 


13.44 


0.0111 ±0.0005 


Ctio/Andicam 


GD 


0.82 (I) 


12.23 


0.033 ± 0.002 


Ctio/Andicam 


GD 


1.23 (J) 


10.44 


0.107 ±0.005 


Ctio/Andicam 


GD 


1.63 (H) 


9.79 


0.124 ±0.006 


Ctio/Andicam 


GD 


2.19 (K) 


9.32 


0.125 ±0.006 


Ctio/Andicam 


GD 


3.60 


8.38 


0.125 ±0.003 


Spitzer/Irac 


M 


4.50 


7.91 


0.123 ±0.001 


Spitzer/Irac 


M 


5.80 


7.42 


0.124 ±0.003 


Spitzer/Irac 


M 


8.00 


6.51 


0.162 ±0.001 


Spitzer/Irac 


M 


24.0 


3.52 


0.280 ± 0.003 


Spitzer/Mips 


S 


7.6-37 




Spitzer/Irs low resolution spectrum 


B,S 


70.0 




0.18 ±0.02 


Herschel/PacsPhot 


BR 


160.0 




0.069 ± 0.007 


Herschel/PacsPhot 


BR 


870.0 




< 0.036 


Apex/Laboca 


NP 


unused data ... 


0.45 (B) 


15.07 


0.00399 


MSSS0 2.3m 


L,BR 


0.558 (V) 


13.97 


0.00940 


SAAO lm 


La,BR 


0.695 (R) 


12.98 


0.0199 


SAAO lm 


La,BR 


0.90 (I) 


11.77 


0.0494 


SAAO lm 


La.BR 


1.24 (J) 


10.51 


0.102 ±0.003 


2Mass 


T,NP 


1.65 (H) 


9.83 


0.125 ±0.004 


2Mass 


T,NP 


2.17 (K) 


9.43 


0.114 ±0.004 


2Mass 


T,NP 


3.80 (U) 


8.14 


0.14 ±0.006 


Vlt/Icsaac 


H,IP 


25.0 




0.298 ± 0.03 


Iras 


I 


60.0 




0.281 ±0.04 


Iras 


I 


70.0 




0.137 ±0.008 


Spitzer/Mips 


S 


160.0 




< 0.12 


Spitzer/Mips 


G 


Detections are 


listed as 


F v + a whereas 


non-detections are 


listed as 



< 3o\ Abbreviations for references, data reduction and flux conversion 
are: GH = Greg Herczeg, this paper; GD = G. Duchene, this paper; BR 
= B. Riaz, this paper; NP = N . Phillips, t his paper; IP = I. Pascucc i, 
this paper; M =lMegeath et al. ( 2005); S = Sicilia ^Aguilar et aD d2 009): 
B =lBouwman et alj j2006h; L^Ilvq et al.l j2004l) ; La = lLawson et alj 
d2002h : H =lHaisch et alJ 120051): T= ICutri et alj (2MASS Point Sourc e 
Ca taloguel2003l); NMoshir et alj (IRAS Faint Source Catalogue 1990); 
G ^Gautieretd1 l l2008l) . 



fluxes adopted in our analysis. The new photometric fluxes are 
substantially lower than previously published fluxes, see Table[2] 

3.3. APEX/LABOCA photometry 

An upper limit at 870 /im has been obtained from new contin- 
uum maps of the rj Chamaeleontis a ssociation taken wi th the 
LABOCA bolometer array on Apex (ISiringo et all |2009). The 
data was reduced using the Bolometer array data Analysis (Bo A) 
software, with a pipeline optimised for faint compact sources. 
Fluxes were extracted within BoA by fitting the amplitude of a 
beam-sized Gaussian at specified positions within a map with 
a pixel scale of 4.6"/px. For each target the flux was extracted 
at points in a 5 x 5 rectangular grid centred on the target with 
a spacing of 37" (twice the beam FWHM). The sample stan- 
dard deviation of the 24 off-source measurements is the 1-cr 
uncertainty quoted here. Using aperture photometry instead of 
beam fitting, with a variety of aperture sizes, and aperture cor- 
rections computed from the beam, yields very similar results 
(9+11 mJy/beam for apertures with r = HWHM). 




0.1 l.o 

A, [um] 



Fig. 2. Observed UV spectra and the stellar atmosphere model 
spectrum for ETCha (see Sect. 14. U . Blue: Fuse spectrum of 
TW Hya divided by 20. Green: Hst/Cos spectrum of ET Cha. 
Orange and red: G230L and G430L Hst/Stis spectra of ET Cha. 
Black: stellar atmosphere model spectrum with T e ff = 3400 K, 
log g — 3.89 and R+ = 0.84 R , at a distance of 97 pc, plotted as 
dotted line (unused) for A <450 nm. 

Table 3. Band integrated UV-fluxes for ETCha. 



spectral band [A] 


integrated spectral flux [erg/cm 2 /s] 


912-1110 


8.03 x 10" 14 


1110- 1450 


1.39 x 10~ 12 


1450 - 2050 


9.73 x 10" 13 



3.4. UV observations 

ETCha was observed in the far-ultraviolet with the Cosmic 
Origins Spectrograph (COS) and in the near-ultraviolet and blue 
with the Space Telescope Imaging Spectrograph (STIS) on board 
of the Hubble Space Telescope (HST) as part of the programme 
"Disks, Accretion and Outflows of TTau Stars" (P.I. G. Herczeg) 
on 5 Feb. 2010. The far-ultraviolet observations were obtained 
with the G130M and G160M gratings with 3891s and 4501s in- 
tegration times, respectively, to cover the 1 140-1790 A spectral 
range with a resolution of ~18 000. The data were reduced with 
the COS calibration pipeline CALCOS and individual segments 
we re combined w i th the IDL coaddition procedure described 
bv iDanforth et alJ (120101) . The near-ultraviolet and blue obser- 
vations were obtained with the G230L and G430L gratings with 
integrations of 45s and 680s, respectively, to cover the 1700- 
5700 A spectral range with a resolution of ~2000. The data were 
reduced with custom-built software in IDL. 

A Fuse spectrum of ET Cha covers the spectral region short- 
ward of 1140 A but has poor signal-to-noise. Instead, the 912— 
1 140 A fl ux was estimated from a high- quality Fuse spectrum of 
TW Hya dJohns-Krull & Herczegll2007h . scaled to the flux in the 
1250-1700 A bandpass. All different data have been smoothed to 
a common low resolution, see Fig. [2] Integrated fluxes over three 
different UV bands as listed in Table [3] 

3.5. Re-analysis of high-resolution optical spectroscopy 

A high-resolution optical spectrum of ETCha was acquired 
on June 26, 2002 with the 3.9m Anglo-Australian Telescope 
Aat and University College London coude echelle spectrograph 
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Table 4. Measured optical emission line properties. 
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Fig. 3. Optical line spectra of ET Cha. All line profiles are con- 
tinuum subtracted and normalised to line peak emission = 1 . The 
forbidden oxygen line [OI] 6300 A is the strongest emission line 
by a factor 5 or more, see Table|4] hence the increased noise level 
in the other lines. The dip in the [OI] 6300 A line coincides with 
the position of the telluric component which has been carefully 
removed, see text. The Lil 6708 A absorption line has been used 
to determine the stellar velocity. 



(Ucles) as first published by lLawson et all d2004l) . A 1.5"slit 
was used delivering a resolving power of ~30000 (lOkm/s at 
6300 A) and covering the wavelength range between 4980 and 
9220 A. The spectrum was calibrated using dome-flats, bias 
frames and ThAr arc frames, making us e of standard library r ou- 
tines such as ccdproc within IRAF (see lLawson et aUl2004 for 
a detailed description). We have paid particular attention to the 
removal of telluric contributions to these lines. We removed the 
OI6300A telluric contribution using the RECX 10 spectrum, a 
star of similar spectral t ype but without any evidences of an ac- 
cretion disk or outflow (iLvo et aUl2003l) . meaning that the mea- 
sured OI6300 A emission from RECX 10 is thoroughly due to 
the atmosphere of the Earth. 

The average air-masses during observations of ET Cha and 
RECX 10 are 2.6 and 2.7, respectively. We made a [OI] 6300 A- 
map by first removing the [OI] feature from the RECX 10 spec- 
trum, and then subtracting this edited spectrum from the orig- 
inal RECX 10 spectrum. The ETCha spectrum was divided by 
the [OI] 6300 A-map to remove the OI telluric line. We used the 
photospheric Lil absorption lines at 6707.76 A and 6707.91 A 
to measure a stellar radial velocity of ~22km/s (~34.6km/s af- 
ter heliocentric correction). 

Figure [3] shows the observed profiles of ETCha in the oxy- 
gen line [OI] 6300 A, with 3 other optical forbidden emission 
lines overplotted that trace outflows, plus a Li I absorption line 
to determine the systematic stellar velocity. The [OI] 6300 A line 
shows a broad component centred around the stellar system- 
atic velocity (low-velocity component LVC), and a blue compo- 
nent shifted by about 42 km/s (high- velocity component HVC). 
We have fitted the HVC and LVC by two Gaussian profiles. 
Measured equivalent widths are listed in Table |4] Equivalent 
widths for these lines ar e converted to line lu minosities using 
the procedure outlined in Hartigan et al.l (119951) . assuming a dis- 



line 


remarks 


EW [A] 


L llnc [10- 7 L o ] 




FWHM 


[OI] 6300 A 


HVC 


-6.0 


170 - 260 


-42 


47 ±15 


[OI] 6300 A 


LVC 


-5.4 


150 - 230 





38+15 


[SII] 6731 A 


all 


-1.8 


50-76 


-35 


n.a. 


[SII] 6716 A 


all 


-0.5 


14-21 


-32 


n.a. 


[Nil] 6583 A 


all 


-0.05 


1.4-2.1 


-40 


n.a. 



EW = equivalent width, HVC = high-velocity component, LVC = low- 
velocity component. The luminosity interval corresponds to the uncer- 
tainty in the red continuum as derived from photometry R = (13.44- 
12.98) mag, compare Table [2] Negative values for the centre velocity 
v cen indicate a blue-shift, n.a. means no values derived. v ccn and FWHM 
are in [km/s]. 
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Fig. 4. v rot sin/ determination: line profile Fourier transforms. 
The Fourier transform of a photospheric line at 6709 A (black 
line) and a telluric line at 9196 A (blue line) are shown. The 
location of the first minimum of the photospheric line profile 
provides a measurement of v lot sin i. The method requires that 
the main line broadening mechanism is rotation velocity, that 
is, the first minimum of the photospheric line must be located 
at a lower frequency than telluric or arc lines, as shown in the 
figure. Arbitrary normalisation constants have been applied to 
the curves to improve readability. 



tance of 97 pc and zero visual extinction. As expected for out- 
flows, the [Nil] 6583 A, [SII]6716A and [SII] 6731 A emission 
lines emphasise the HVC, and are all blue-shifted by about the 
same margin. 

3.6. Projected stellar rotational velocity 

The high resolution spectrum presented in Sect. 13.51 was addi- 
tionally used to determine the projected rotational velocity of 
the star v lot sin i. This quantity is typically derived from empir- 
ical relations between the full-width at half-maximum of stel- 
lar ab sorption lines and v rot sin i (see e.g. Martin ez- Arnaiz et all 
2010). How ever, in this paper we use the Fourier transform of the 
line profile dGrav. D. F.LU992T) . The full power of this method is 
reveal ed in case of very high s pectral resolution and signal to 
noise (Re iners & Schmitn, 120031) . but has als o been applied suc - 
cessfully to moderate and low quality data dMora et all 1200 lb . 
using the following simplifications. 
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Fig. 5. Bayesian probability distributions for effective temperature T e g, visual extinction Ay, distance D and stellar radius /?+ for 
ETCha using our new simultaneous BVRIJ photometry (solid histograms). The vertical dotted lines bracket the range of values 
explored for each parameter while the vertical dashed line represent the absolute best fit photospheric model. 



The location of the first minimum of the Fourier transform 
of any isolated photospheric line profile provides a measurement 
of v sin i. If many lines are available, the average value and stan- 
dard deviation can be estimated. The advantage of this method 
is that the measurements are independent of empirical calibra- 
tions and synthetic spectra, provided that enough lines can be 
identified and the dominant line broadening mechanism is stel- 
lar rotation. In terms of Fourier analysis, the latter requirement 
is equivalent to the following condition: the first minimum of 
the Fourier transform of photospheric line profiles must be at a 
smaller frequency than the Fourier transform of lines not affected 
by rotation velocity (e. g. arc or telluric lines). In addition, the al- 



gorithm is robust against the introduction of noise, and capable 
of working with low signal to noise data. 

The method has been applied to the visible Ucles spectrum 
of ET Cha as described in Sect. 13.51 Few lines could be used 
for the analysis, due to the low signal to noise ratio (~10 for 
unbinned data). A synthetic Kurucz spectrum with physical pa- 
rameters similar to those listed in Sect. I4.1l was used to identify 
isolated atomic photospheric lines. After clipping bad lines, only 
16 lines remained in the final average used below. Figure|4]shows 
a sample Fourier transform of both a photospheric and a telluric 
line to illustrate the method. We obtain 




^4oO -50 50 100 

Vlsr (km/s) 

Fig. 6. Deep APEX spectrum of ETCha at 345.796GHz 
(866.963/mi), showing the mean beam temperature T m \, versus 
local standard of rest velocity vi sr . The 12 CO J = 3 — > 2 line was 
not detected. 



v rot sin i = (15.9 ±2.0) km/s. (1) 



3.7. APEX n CO 7 = 3— > 2 

The source was studied in the J — 3 — > 2 rotational transition of 
I2 CO at 345.796GHz using the Apex telescope and heterodyne 
receiver system twice, in April and July 2010, for a total on- 
source time of 100 minutes. A standard on-off switching scheme 
was used, with reference position 180" off in Azimuth. With an 
aperture efficiency of j] m b = 0.6, the effective area of the tele- 
scope is 67.9m 2 , and the beamsize Q.a at 345.796GHz is 13.7". 
The CO J — 3 — > 2 line is commonly detected from disks of radii 
100 AU or more; however, in the case of ETCha, the line was 
not detected (Fig|6]l, down to an rms of 17 mK (antenna temper- 
ature) for 0.425 km/s wide velocity bins. Assuming the under- 
lying linewidth of the emitting region is 7.5 km/s, as predicted 
by the best fitting model (Sect. [5} and typical of CO lines from 
disks, the upper limit to the line flux is 3cr = 5.5 x 10~ 20 W/m 2 . 
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4. Disk modelling 

The strong near-mid IR continuum excess of ET Cha (see Fig. [7} 
as well as the narrow H2 v— 1 — > 0S (1) emission line at 2.122/im 
( Rams ay Howat & GreavesL 120071) indicate the presence of a 
gas-rich protoplanetary disk. The optical high-resolution spec- 
trum with blue-shifted emission lines shows, however, that in 
addition to the disk, there also exists an outflow that might con- 
tribute to the gas emission lines. In the following, we present 
a combined gas+dust disk model to explore in how far a disk 
model alone can explain all photometric and spectroscopic ob- 
servations of ET Cha. The possible contribution of an outflow is 
discussed in Appendix l6.4l 

The disk modelling procedure in this paper consists of three 
phases. In phase 1, we make a Bayesian analysis, based on the 
new BVRIJ photometric data, to obtain best fitting values for 
the stellar parameters. In phase 2, we fix a few more parameters 
like the inner disk radius based on physical assumptions, e.g., 
according to stellar luminosity and dust sublimation tempera- 
ture, and in phase 3, we use the thermo- chemical code ProDiMo 
dWoitke et ali r2009: Ka mp et all l2009) to calculate the gas- and 
dust temperature structure in the disk, the chemical and ice com- 
position, and to fit all remaining disk, dust, and gas parameters 
of the model to the observations as good as possible. 

4.1. Stellar parameters 

While our new JHK photometry is consistent with the 2 Mass 
photometry, the VRI photo metry is about 0.5 - 0.7 mag fainter 
than the data published in dLawson et all l2002h . Indeed, our si- 
multaneous photometry appears to represent the 'faint' state for 
ET Cha which, in the accretion hotspot scenario, is purely pho- 
to spheric. 

To estimate the stellar properties, we therefore adopt our new 
simultaneous photometry. We only fit the BVRIJ fluxes to en- 
sure that our estimates are not biased by the near-IR excess from 
the disk. We computed a grid of photospheric mo dels, varying 
r e ff ( using the logg = 4.5 NextGen m odels from (Allard et all 
1997), A v (using the R v = 3. 1-law from lCardelli et all 1 19891) . ~kZ 
and the distance d to ETCha within generous ranges. The opti- 
cal extinction is defined by Ay = 2.5 log 10 (e) Ty where Ty is the 
optical depth in the visual due to interstellar dust extinction. We 
note that the last two parameters are degenerate. 

The results are shown in Fig. [5] A very good fit is obtained 
with J eff = 3300 - 3400 K, d = 91 - 100 pc, R* = 0.84-0.89 R @ 
(corresponding to L+ = 0.085 L ) and Ay < 0.2 mag. This is 
in exc ellent agreement with the findings of ILuhman & Steeghsl 
(2004), including the effective temperature which they derived 
exclusively from their spectrum and not from photometry. For 
the disk modelling, we adopt T eff = 3400 K and R+ = 0.84 R G , 
which is both the most probable combination of parameters and 
the closest to the spectroscopically derived stellar properties. 

Finally, we note that conducting the same analysis on the 
previous photometric dataset (obtained in the "bright" state) 
yields both a hotter central star (r e ff a4000K) and a substantial 
foreground extinction (Ay = 0.75 mag), both of which are incon- 
sistent with our prior knowledge of the source (see Sect. |2}- 

In the following, we use the nearest NextGen stellar in- 
put spectrum for T e ff = 3400 K, L+ = 0.085 L and solar 
metallicities (resulting in log g = 3.89 and /?* = 0.84 R Q for 
M+ — 0.2 Mo). A ccording to the stellar evolutionary models of 
Siess et al. (2000), our choice of T e g and is consistent with a 
6-8 Myr old star of mass (0.2-0.3) M and solar metallicities. 



4.2. Stellar UV-excess 

In the UV, stellar activity and mass accretion create an ex- 
cess emission with respect to classical stellar atmosphere mod- 
els. Therefore we use our composite observed UV spectrum of 
ETCha (Fig. [2j at /l<450nm as input instead. Although strong 
Lya emission from ET Cha was detected with Hst/Cos, no emis- 
sion was detected near line centre because of neutral hydrogen 
in our line of sight to the star. Based on an analysis of a sim - 
ilar Lya emission profile from TWHya dHerczeg et all |2004), 
we roughly estimate that the total flux in this line is twice the 
detected value. Therefore, we have doubled the detected flux be- 
tween 1210 A and 1220 A for the construction of our UV input 
spectrum "as seen by the disk". This yields an integrated flux of 
2.22 x 10~ 12 erg/cm 2 /s in the 1 1 10-1450 A band, in contrast to 
Table [3] We have converted these integrated UV fluxes into ad- 
ditional photometry points, see Table [2] and plotted these points 
with large errorbars in Fig. [7] These data result in a fractional 
UV luminosity /uv = Luv/L* - 0.017 with the UV luminos- 
ity Lpy being integr ated from 9 1 2 A to 2500 A as introduced by 
Woit ke et alT(l2010h . The "extra UV" is not important for the 
spectral energy distribution (SED) modelling, but is essential for 
the gas modelling, because the UV irradiation is a decisive factor 
for the gas heating and chemistry in the disk. 

4.3. Disk inclination and inner radius 

For the disk modelling in Sect. 14.51 we have fixed two parame- 
ters to reduce the dimension of parameter space for the fitting 
problem, namely the disk inclination and inner radius. 

First, the disk inclination is assumed to be i = 40° as mea- 
sured from face-on orientation. As long as the disk is not inter- 
secting the line of sight to the star, this parameter does not have a 
large impact on the model results concerning both SED-analysis 
and line flux predictions, and we do not have clear observations, 
e.g., images, that would allow us to determine this quantity un- 
ambiguously. Large values fo r the disk inclination are supported 
by the H ff line analysis (60°, Law son et all 120041) . the low out- 
flow velocities observed (Sect. 13.5b , and by the strongly rotation 
broadened stellar absorption lines (Sect. 13.6) , favouring an edge- 
on rather than a face-on disk orientation. 

However, our disk modelling suggests large disk scale 
heights in the inner disk regions, needed to reproduce the strong 
near-IR excess (see Fig. lC.41 >. For inclinations of 60°and higher, 
the observer's line of sight to the central star would therefore in- 
tercept the disk. This would cause a dramatic reduction of the ob- 
served fluxes at optical and UV wavelengths, as well as a strong 
reddening of the optical colours. This is inconsistent with our 
initial assumption that the star is seen directly, used to estimate 
the stellar properties. In principle, one could increase both the 
intrinsic stellar luminosity and effective temperature of the star 
to compensate exactly for this reddening, but to match the ob- 
served SED, the required T e ff is well beyond values consistent 
with the M3-M3.5 spectroscopic classification of ETCha. For 
inclinations ^ 60° we furthermore observe that the 10 /im and 
20 fim silicate features are no longer in emission, which clearly 
contradicts of observed SED of ET Cha. 

Second, concerning the disk inner radius, the strong and con- 
tinuous near-IR fluxes of ET Cha suggest a disk which extends 
inward close to the star, with no dust holes or gaps. T Tauri 
disks are generally assumed to be truncated by the stellar mag- 
netosphere near the co-rotation radius, with material accreting 
along magnetic field lines ont o high-latitude regions of the star 
(lKoniglL[l99TI:IShu et allll994h . From the analysis of optical ab- 
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Fig. 7. Comparison of model SED to observations. Blue bullets 
are photometric data points used for the fitting (see Table |2). 
Light blue symbols refer to other, unused photometry points. 
Symbols with downward arrows indicate 3<x upper limits. The 
black line is the SED from our best-fitting model. The red line 
shows the assumed stellar input spectrum, completed by piece- 
wise constant fluxes F v in three UV bands as obtained by inte- 
grating over our UV observations (Table [3}. 



sorption lines of ET Cha (Sect. 13.6b we have inferred a projected 
stellar rotation velocity of v rot sin i = (15.9 ±2.0) km/s. At / = 40°, 
this translates into v rot = (25 ± 3) km/s (about 12% of the break- 
up velocity) and, by assuming M* = 0.2 M and R+ = 0.84 R e 
(Sect. l4~Tl ). into a stellar rotation period of P- (1.7 + 0.2) days, 
and a co-rotation radius of (0.015 + 0.002) AU. 

However, the dust temperature at the co-rotation radius 
(0.015 AU) is larger than 1500 K, associated with the sublima- 
tion temperature of silicates. Therefore, we have assumed a 
slightly larger inner disk radius, R\ n = 0.022 AU, where the dust 
temperature is about 1500 K. 

4.4. A single disk model 

For a single disk mo del, we use the radiation thermo -chemical 
disk code ProDiMo (IWoitke et all 120091: iKamp et all l2009h to 
calculate the dust continuum radiative transfer, and the gas ther- 
mal balance and chemistry throughout the disk. We use 10 ele- 
ments, 76 gas phase and solid ice species, and 992 reactions in- 
cludin g a detailed treatment of UV-photorates (see Kam p et all 
2009), H2 formation of grain surfaces, vibrationally excited H? 
chemistry, and ice formation (adsorption, thermal desorption, 
photo-desorption, and cosmic-ray desorptio n) for a limited num- 
ber of ice species (see Woitke et al., 2009, for details). We also 
use ProDiMo in this paper to compute all observables includ- 
ing the SED, images, and gas emission line fluxes and pro- 
files. Latest improvements to the ProDiMo model include X- 
rays chemistry and heating, a parametric prescription for dust 
settling, UV fluorescent pumping, PAH ionisation and heat- 
ing/cooling, [OI] 6300 A pumping by OH-photodissociation, 
H2-pumping by its formation on grain surfaces, formal solutions 
of the line transfer problem, and chemical heating. These im- 
provements are explained in the Appendices IA.2l to lA. 81 



Our modelling of ET Cha is based on a prescribed disk den- 
sity structure, using power-laws for the surface density distribu- 
tion and disk flaring, see Appendix lA.il This approach allows for 
a rapid mod el computation (avoi ding the structure iteration loop, 
see Fig. 1 in Woitk e et alll2009l) . and is hence more suitable for 
a deep search for fitting values in parameter space. In this mode, 
the disk code uses altogether 25 physical parameters, most of 
which are considered as fixed for the modelling of ET Cha (for 
instance the stellar properties, see Table |5}- 

4.5. Parameter fitting procedure 

The determination of the remaining free disk, gas and dust pa- 
rameters, like the total disk gas mass for instance, is a key ob- 
jective of the modelling of individual targets. These parameters 
are "determined" in this paper by considering the values (better 
the ranges of values) that lead to a satisfying match between pre- 
dictions and observations in the frame of the model, henceforth 
called "solutions". Practically all determinations of properties of 
astrophysical objects work that way, even if it's just a simple 
one-dimensional dependence between property and observable, 
because the desired quantities are rarely accessible via direct ob- 
servations. Thi s is the wel l-known problem of model inversion 
in astronomy with all the usual shortcomings and 

concerns which can be subdivided into four families: 

1. Concerns about the quality of the model itself (missing 
physics, poorly determined input quantities like cross- 
sections and rate coefficients, numerical issues like grid 
resolution, etc.). 

2. Concerns about the completeness of solutions found in 
parameter space, for instance, how can we be sure of having 
found the best solution? 

3. Under-determination, i.e. a weak dependence of observ- 
ables on model parameters in combination with considerable 
uncertainties in the observables. 

4. Non-uniqueness of inversion procedure (multiple solutions) 

Given the complexity of the disk model at hand, with 8-14 free 
parameters for the modelling of the disk of ETCha, we have 
to take into account the possibility that different solutions exist 
which practically achieve the same degree of agreement between 
model results and observations. Whereas in a low-dimension pa- 
rameter space this seems to a be a weird, seldom special case, it 
occurs frequently in high dimensions, with numerous local min- 
ima. The manifold of solutions in a multi-dimensional parameter 
space can be (and usually is) astonishingly rich and complicated 
in structure. 

Furthermore, it is important to note that, in a high-dimen- 
sional parameter space, an exhaustive search is practically im- 
possible. This is even more so for our disk modelling, since one 
complete model takes about 1 CPU hour on a single 3 GHz pro- 
cessor machine. For 10 parameters, with 20 well-selected values 
around a main solution each, one would need to run 20 10 « 10 13 
models which would take about 1 .2 x 10 9 CPU yrs. 

Therefore, we are not able to fully devitalise the concerns 
(1) to (4), but have to face the fact the any parameter determi- 
nation by model inversion in high-dimensional parameter space 
is an intrinsically uncertain business. Our strategy in this paper 
is as follows. We use an evolutionary strategy to find a hand- 
ful of well-fitting disk models, corresponding to different local 
minima in parameter space. Among these, we select a "best- 
fitting" model, the properties of which are described in detail in 
Appendix ICl The best-fitting model is then re-run with different 
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choices of input physics in Appendix[D] to explore the principal 
effects of e. g. X-rays and the treatment of ^-formation. 

We then conclude about the confidence intervals of derived 
parameter values in Sect.|5]by considering (i) small deviations of 
single parameters around the best fitting solution, (ii) overall ex- 
perience from fitting by hand and variance of different solutions 
found by the evolutionary strategy, and (iii) direct constraints 
from physical arguments like optical depth and column densi- 
ties. We admit that this modelling procedure is not entirely sat- 
isfactory. A more complete discussion of our modelling strategy 
will be the topic of a forthcoming paper. 

4.6. Fit quality and evolutionary strategy 

Mathematically, all model inversion techniques can be formu- 
lated in terms of a certain strategy to minimise x 2 , i- e. to find a 
minimum of the deviations between model predictions and ob- 
servations in parameter space. In this paper, we use the following 
logarithmic measure of these deviations as 



Table 5. Parameters of the best-fitting disk model. 



with 



^ mod 



3 a* 



(2) 



F' ohs and cr' F are an observed flux and its uncertainty, and F' d is 
the predicted flux by the model. The logarithmic nature of our^- 2 
is motivated by the need to assign an equally large (bad) number 
to x if a model flux is a factor 10 too large as if it is a factor 10 
too small, analog to considering deviations in magnitude^. 

Equation (f2| is applied separately to all photometric data, all 
/l-points in the Spitzer spectrum, and to the line fluxes of CO J = 
3^2, [OI]63^m, [OTJ6300A (LVC) and o-H 2 2. 122 /mi, re- 
sulting in*p hot , xj pit , and^- 2 ine , respectively. The total fit quality 
of a model is then calculated as 



X 1 = w Pho ^ ho 



pit 



+ w L ™xl 



(3) 



F Phot , w Splt and w Llne , normalised as 



with adjustable weights w 

VyPhot +>v Spit + ^Line^ 

We have applied th e (1,11)- evolutiona ry strategy with adap- 
tive step-size control o flRechenber g (2000) to minimise^, i. e. to 
find best-fitting values of our remaining free gas, dust and disk 
parameters. The strategy uses 1 parent producing 11 offsprings 
with slightly modified parameters, the best of which will become 
the parent of the next generation (the parent always dies). The 
step-size 6 is transmitted to the children and treated as addi- 
tional parameter to be optimised. After some experiments, we 
have chosen w phot = 0.35, w Spit = 0.45 and w Line = 0.2 for opti- 
mum performance of the evolutionary strategy. 

Figure [8] visualises one exemplary run of the evolutionary 
strategy, showing the changing model parameters and results 
over 300 generations. We started from a generic disk setup with 
M gas = 1O~ 3 M , gas/dust ratio of 100, and an outer disk radius of 
^out = 100 AU. This model fails badly to explain the SED and, 
in particular, the line fluxes. CO J — 3 — > 2 is about 10 times too 



Defining a magnitude as m' obs oc logF' obs , the error thereof is cr' m 



quantity 


symbol 


value 


stellar mass 


M+ 


0.2 M Q 


effective temperature 


T eff 


3400 K 


stellar luminosity 


U 


0.085 L 


disk gas mass* 




6.1 x 10~ 4 M o 


inner disk radius 




0.022 AU 


outer disk radius* 


^out 


8.2 AU 


column density power index* 


£ 


-0.020 


reference scale height* 




0.011 AU 


reference radius 




0.1 AU 


flaring power index* 


P 


1.09 


disk dust mass* 




2.6 x 10~ 8 M 


minimum dust particle radius 




0.05 yum 


maximum dust particle radius 


n 

"max 


1 mm 


Hnst si7P H i ct nnwpr lnrlpx 


y 


4.1 


minimum settling particle size 


d 


o 


dust settling power index 


s 





dust material mass density 


Pgr 


3 g cm -3 


dust composition 


Mg 2 Si0 4 


32.9% 


(volume fractions) 


amorph. carbon 


24.4% 




MgFeSi0 4 


23.0% 




SiO, 


8.8% 




MgSi0 3 


7.6% 




cryst. silicate 


3.3% 


strength of incident ISM UV 


ISM 

A. 


1 


cosmic ray H 2 ionisation rate 




5 x 10~ 17 &- [ 


PAH abundance rel. to ISM* 


/pah 


0.081 


chemical heating efficiency* 


yChem 


0.55 


a viscosity parameter 


a 





disk inclination 


i 


40° 


distance 


d 


97 pc 



fllogCF' )/dF' 



o~' F /F' h . and hence the deviation A, = (m'. 



m' is equivalent to Eq. l|2j. 



Parameters marked with * have been varied in the evolutionary optimi- 
sation run depicted in Fig. [8] The values of the other parameters have 
been assumed, fitted by hand, or have been obtained from additional 
evolutionary optimisation runs not discussed here. 



strong and both [OI] 6300A and o-H2 2.122/mi are more than a 
factor of 100 too weak. However, after about 300 generations, 
the model has achieved a good fit of all available continuum 
and line observations. The final parameter values are far from 
these initial guesses, yielding a disk that is less than 10AU in 
radius and much less massive. All models from about gener- 
ation 50 onwards fit the dust observations about equally well, 
the SEDs from these models are almost indistinguishable by eye 
when plotted as in Fig. [7] A good fit to the line observations was 
achieved only from generation 200 onward, after the disk radius 
has shrunk to about R oa $10 AU while continuing to fit the dust 
observations about equally well. Thus, the line observations can 
help to break the degeneracy of SED-fitting. 



5. Results 

We identify the result of the evolutionary run depicted in Fig. [8] 
as our main "best-fitting" model. The resulting parameters are 
listed in Table [5] the computed line fluxes are summarised in 
Tableland the resulting SED is plotted in Fig. [7] More details 
about the internal physical and chemical structure of this disk 
model are shown in AppendixICl 

However, the identification of a best-fitting model was not 
at all straightforward. Altogether 17 runs of the evolutionary 
strategy have been executed with about 50 to 300 generations 
each, choosing different parameters to be varied, different ini- 
tial guesses of the model parameters, or using different setups 
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Fig. 8. Auto-optimisation of model parameters to fit all continuum and line observations of ET Cha. The first row shows the evo- 
lution of some model parameters as function of generation (total gas mass M gas [M ], total dust mass Md us t [M Q ], outer disk radius 
^out [AU]). The two l.h.s. plots in the second row show the resulting x f° r photometry and Spitzer spectrum. The other four plots 
shows the resulting line fluxes as labelled, where the orange bars indicate the observed line flux ranges with 1 <x errors. At generation 
94, the evolutionary strategy got stuck in an unsatisfactory solution (local minimum) with well-fitting SED, but very bad line fits 
A'Line ~ 5, i. e. a mean deviation between model and observed line fluxes of (e 5 »150) cr. We have manually changed parameters 
at generation 95, increasing the disk mass by a factor 10. The evolutionary strategy then found a well-fitting model which can 
reproduce all available line and continuum observations within about 1 cr on average. 



of the evolutionary strategy. Not all of these runs succeeded to 
find well-fitting solutions. Some almost equally well-fitting so- 
lutions are listed in Table|7] Based on these diverse model inver- 
sion results, we have to be very careful with conclusions about 
the disk properties of ETCha. The following section summarises 
our confidence intervals for the various disk shape, dust and gas 
parameters and discusses which observations are key for their 
determination. 



5.1. Dust mass 



ET Cha's spectral energy distribution (SED) is characterised by a 
strong near and mid IR excess relative to the star, similar to other 
much more massive T Tauri disks, see Fig. [7] However, in the 
far-IR, the fluxes are very faint. 69mJy at 160,/im was too weak 
to be detected by Spitzer, and we have not detected the object at 
870 yizm with Apex. Thus, the new Herschel photometry points 
at 70 fim and 160 /mi allow for the most complete SED analysis 
of the source to date. 
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Table 6. Computed line fluxes [1CT 18 W/m 2 ], at d = 97 pc and 
i = 40° from the best-fitting disk model, in comparison to the 
observations. Non-detections are listed as <3cr upper limits. 



line 


/t l//mj 


observed 


model 


[OI] J Pi -> J P 2 


63.18 


30.5 ± 3.2 


34.5 


[OI] To -» 3 Pi 


145.52 


< 6.0 


2.6 


[OI] D 2 — > P2 (HVC) 


0.6300 


73 ± 25 


— 


[OI] 'D 2 -» 3 P 2 (LVC) 


0.6300 


65 ± 25 


69.6 


rrrn 2 p, ,„ — » 2 p, ,„ 


1 ^7 74 


< ■ y A} 


u.l 1 


C07 = 3^2 


866.96 


<0.05 


0.014 


CO/ = 29^28 


90.16 


< 9.6 


4.9 


CO/ = 33^32 


79.36 


< 24 


3.3 


CO/ = 36^35 


72.84 


< 8.0 


2.6 


o-H 2 v=1^0S(l) 


2.122 


2.5 ±0.1 


2.4 


o-H 2 2 2 , — >2j 2 


180.49 


< 5.2 


1.1 


o-H 2 2 12 ^ lot 


179.53 


< 5.0 


1.4 


o-H 2 4 32 ->3 12 


78.74 


< 30 


11.1 


p-H 2 3 2 2~ >2n 


89.99 


< 9.6 


6.4 



Table 7. Parameters of different, about equally well-fitting mod- 
els from different runs of the evolutionary strategy. 



parameter 


model 1 


model 2 


model 3 


model 4 


M gas [10- 4 M o ] 


0.088 


0.65 


6.1 


25 


M dust [10- 8 M G ] 


3.3 


2.6 


2.6 


3.7 


* ou, [AU] 


5.9 


7.5 


8.2 


8.3 


H [AU] 


0.0103 


0.0096 


0.0110 


0.0108 


6 


1.16 


0.046 


-0.020 


0.008 


P 


1.33 


1.15 


1.09 


1.07 


P 


3.2 


3.9 


4.1* 


4.2 


a s [p:m] 


0.05* 


0.01* 


0* 


0* 


s 


0.01 


0.13 


0* 


0* 


/pah 


0.081 


0.098 


0.081 


0.081 


yChem 


0.14 


0.20 


0.55 


0.09 


results 


XPhot 


0.48 


0.78 


0.50 


0.45 


XSpit 


0.91 


0.51 


0.51 


0.50 




1.03 


1.04 


1.25 


1.02 



Each final model is located in a solid local ^--minimum, the essence 
after trying a few thousand models. Parameter symbols are explained in 
Table [5] Other fixed model parameters are as listed in Table [5] Model 3 
was selected as best-fitting model. Parameters marked with * have been 
fixed during the evolution. 



In all computed disk models for ET Cha, the disk is optically 
thin at 160 /mi. The observed flux at distance d is hence given by 

_ K? s B v ({T dust ))M dust 



where (7d us t) is the dust mass averaged dust temperature, the 
dust absorption coefficient in [cm 2 /g(dust)] and B v the Planck 
function. In the best-fitting model we measure (7d us t) ~ 56 K, 
and Ky bs = 34cm 2 /g(dust) at 160 /mi, which, according to 
Eq. ©, results in a 160 /im flux of 48mJy which is in good 
agreement with both, the computed flux from the full disk model 
(5 1 mJy) and the observed value of 68 mJy. 

Therefore, the observed 160 //m-flux leaves no doubt that the 
mass of the dust in the disk of ET Cha must be low, Md us t ~ 
2.6 x 10~ 8 M according to our best fitting model. More massive 
disks produce too strong 160 /im continuum fluxes. Although 
there are certainly some ambiguities in the dust mass determi- 
nation, e. g. if other dust size parameters are used, leading to a 
different Kl hs , but among all parameters, Md us t certainly belongs 
to those which are only little influenced by other parameters and 



can be determined with some confidence. Across all SED-fitting 
disk models, with varying dust composition and size parameters, 
the value for M dust ranges in (2 - 5) x 10~ 8 M Q . 

5.2. Dust disk characteristics and particle sizes 

According to the low dust mass we derive, the disk of ET Cha has 
only small optical depths across the disk. Our best-fitting model 
has vertical dust absorption optical depths <1 for A ^ 4 /mi, ex- 
cept for the optically thick 10/mi and 20 /mi silicate features. 
This is valid for all radii since the column density exponent 
e ss 0. However, the radial optical depths are much larger. The 
midplane radial dust extinction optical depth at 1 /mi is about 
150. Therefore, the disk of ETCha is on the borderline between 
optically thin and thick. It has a complex dust temperature struc- 
ture due to radial shielding effects (see Fig. IC.2t . but the vertical 
dust emission could be treated in the optically thin limit at most 
wavelengths. 

Some conclusions about the dust particle sizes can be drawn 
from the shape of the 10 /mi and 20 pm silicate features, which 
are clearly seen in emission for ET Cha (FigJT}. We need a large 
amplitude of opacity variation across these features (see Fig lC. lb 
to model the SED of ET Cha, which favours small, (sub-) micron 
sized dust particles. Since the peaks are optically thick, warm 
and small grains must be located in front of cooler dust along 
the line of sight. We notice that it is easier to fit the observed 
silicate features with small values of the column density power- 
law index e, i. e. a roughly flat surface density distribution, and 
small disk inclinations. For larger e or larger disk inclinations 
(closer to edge-on), the silicate features weaken and eventually 
vanish. 

If dust and gas are well-mixed (no dust settling, e. g. mod- 
els 3 and 4 in Table [7]), the models clearly favour very small 
particles. We can fit the entire SED with a uniform dust pop- 
ulation that is either truncated at about 1 pm, or, alternatively, 
with a continuous size distribution ranging from 0.05 /mi to 
1 mm, but with an unusually large power-law index of p = 4.1. 
The second approach allows for slightly better SED-fits. These 
findings with ProDiMo have been carefully chec ked against 
the Monte C arlo radiative transfer code MCFOST (tPinte et all 
120061 120091) . showing a very good agreement in calculated dust 
temperatures and continuum fluxes. 

However, if dust settling is taken into account (in the ap- 
proximate way explained in Appendix IA.2l i. Table Q demon- 
strates that smaller p — 3.9 or even p — 3.2 are also possible, 
in which case the volume-integrated dust size distribution in 
ETCha would not be unusual at all, close to the default value 
of p — 3.5. Our honest conclusion about ET Cha is hence that it's 
dust must predominantly be made of (sub-) micron particles at 
the disk surface, where the silicate emission features form. We 
have therefore decided to refrain from quoting any errorbars to 
our results concerning the dust size parameters. 

5.3. Dust composition 

In our best fitting model, we assume a n effective mix o f 
about 33% amorphous fosterite Mg 2 SiO^ dJager et all l2003h . 
24% amorphous ca rbon (|Zubko et all 120041) . 23% amorphous 
olivine M gFeSiQ4 ( Dorschne r et all 1 19951) . 9% amorphous sil- 
ica SiOn (|Posch et all 120031) . 8%~amorphous enstatite MgSi0 3 
(Dorschn er et all [19 95 ). and 3% crystalline fosterite Mg2Si04 
(Servoin & Piriou, 1973). The citations indicate our sources 
for the complex refractory indexes of the various pure materi- 
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Fig. 9. Systematic variation of single parameters around the values of the best-fitting disk model, showing the dependencies on total 
disk gas mass M gas [M ] (l.h.s.), on total dust mass Mdust [M Q ] (middle) and on outer disk radius R out [AU] (r.h.s.)- The full black 
line shows the total x, and the dotted, dashed and dashed-dotted lines show its constituents, jfphot, ^Spit and^jjne, see Eq. (f3]l. The 
computed line flux ratios with respect to the observations F\i ns /F \, s (F\i ne /(3cr) for CO 3 — > 2) are shown by the coloured lines as 
labelled, with tickmarks on the right axis). The grey shaded box marks the parameter-interval where x is n °t larger than twice its 
minimum. 
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Fig. 10. Comparison between predicted SEDs of two models 
with different dust material composition. The black model shows 
once more the SED of our best fitting model, with dust compo- 
sition as listed in Table [5] In the green model, we have omitted 
the amorphous carbon. All model parameters are identical oth- 
erwise. 



als. The dust absorption and scattering opacities ar e calcu lated 
by applying effective mixing theory (Bruggemanl Il935h and 
Mie theory, based on these optical data and volume fractions, 
and assuming that the chemical dust composition and unsettled 
dust size distribution is unique throughout the disk. The inclu- 
sion of crystalline fosterite is motivated by the fine- structure of 
the ob serve d second silicate feature at 2 //m dBouwman et all 
2006), see (Sicilia-Agui lar et all 120091) . which shows several 
na rrow peaks close to the fosterite peak positions in the data 
of Se rvoin &Pirioul (119731) . The resulting volume fractions are 
a by-product of our automated fitting procedure from additional 
runs of the evolutionary strategy not shown in Table Q We do 
not claim, however, to have determined the dust composition of 



ET Cha, as our method is focused on fitting the overall shape of 
SED rather than individual dust features. 

The inclusion of 24% amorphous carbon, however, was an 
important step to understand the SED of ET Cha. A comparison 
model without amorphous carbon (green line in Fig. [10} demon- 
strates its impact on the SED. First, amorphous carbon reduces 
the dust albedo at UV to near-IR wavelengths (see Fig. IC.U . 
Pure laboratory silicates have an albedo of about 80% -99% 
around 1 yum. This leads to a substantial starlight amplification 
via scattering by the disk at UV to near-IR wavelengths, which is 
inconsistent with the photometric data. Second, amorphous car- 
bon enhances the absorption and thermal emission in the near- 
mid IR. As a consequence, the disk is more effectively heated by 
star light, and produces more thermal emission shortward of the 
10/im silicate feature, just where ETCha is very bright, result- 
ing in a much better fit of the SED if we include amorphous car- 
bon. However, any other kind of impurities or inclusions, metal- 
lic iron for instance, would cause a similar increase of the dust 
absorption opacities at optical to near-IR wavelengths, amor- 
phous carbon is just one of the options. The formation of small 
alien inclusions in the du st m aterial has be e n de monstrated by 
IHelling & Woltkel d2006l) and IHelling etail 420061) to be a nat- 
ural consequence of the refractory dust formation process in 
somewhat different oxygen-rich environments, namely in brown 
dwarf atmospheres. Concerning ETCha we conclude that strong 
near-IR dust absorption opacities are needed to fit the SED, sub- 
stantially stronger than those of pure amorphous or crystalline 
silicates which have a "glassy" character. 

5.4. Outer disk radius 

Any outer disk radius, at least up to 200 AU, works fine to fit 
the SED alone. However, the non-detection of CO J = 3 — > 2 
and the modest [OI] 63 fim line flux put severe constraints on 
R ont . As Figs. [8] and [9] demonstrate, the CO line flux depends 
very strongly and robustly on R out . All calculated models with 
R out > 25 AU would violate the 3<x upper limit of CO J = 3 — > 2. 
Only models with R out < 10 AU are consistent with the l<x CO 
line flux upper limit. The CO lines are extremely optically thick, 
T"Mne ~ 10 5 (see Sect. 15.5b . Therefore, even if the CO abundance 
was reduced by a factor of 10 3 throughout the disk, for example 
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Fig. 11. Line profiles, intensity maps, and spatial origin of some emission lines according to the best-fitting disk model. The upper 
row of figures shows the calculated continuum level, line profile, and continuum-subtracted frequency-integrated line intensity- 
map for [OI] 63.2 yum (left), [OI] 6300 A (middle) and o-H2V=l— >0S(l)2.122/um (right), based on formal line transfer solutions at 
inclination = 40°. Note the different scaling of velocity axis for the different lines. The lower row of figures visualises the vertical 
line and continuum optical depths as function of radius, the cumulative line flux, and the species density [cm -3 ]. The vertical dashed 
lines indicate where the line flux is coming from, bracketing 70% of the cummulative line flux in radius. The thick black lines in the 
lowest boxes mark those grid cells that contribute most to the vertical line flux; they mark the geometrical depth of the line formation 
region. For the analysis figures in the lower row, we have directly applied the upwards escape probability formalism (without formal 
solutions of line transfer). These analyses are hence only approximate in nature and valid for strictly upward line propagation only 
(inclination = 0°). The line fluxes calculated in this way deviate from the proper results under 40° by less than 8%. 



by very efficient CO ice formation, the conclusion that the disk 
of ET Cha must be small would still be valid. The relatively ro- 
bust determination of R out *(6 - 9) AU is demonstrated by Fig. [9] 
and Table|7] Smaller values of R out are actually inconsistent with 
the measured [OI] 63 yum line flux (unless interpreted as origi- 
nating form an outflow, see Sect. 16.4b . because the outer disk is 
responsible for the [OI] 63 fim line emission (see Sect. 15.5b . 

5.5. Spatial origin and characteristics of gas emission lines 

Before we continue to conclude about the determination of the 
disk gas mass of ET Cha, we first have to clarify where the ob- 
served spectral lines come from and what they tell us. Figure [TT] 
shows the calculated line fluxes and profiles of the three detected 
lines [OI]63jum, [OI]6300A (LVC) and o-H 2 2.122yum. We 
also plot here the vertical line and continuum optical depths and 
the cumulative line flux as function of radius, to facilitate the 
discussion of the lines' spatial origin. These plots are based on 
the best-fitting disk model, but the drawn conclusions about the 



spatial origin of the spectral lines are quite general and valid 
for all calculated disk models that fit the observations of ET Cha. 

The far-IR [OI] 63fim line probes the outer disk layers, about 
2-7 AU in this tiny model, and is optically thick with vertical 
optical depths rii ne w 100-1000 in this radial region^. The con- 
tinuum is optically thin. Since the line is collisionally excited, 
with an excitation energy of about 228 K, the line flux probes 
first and foremost the existence of warm gas (r gas ^ 50 K) in 
the disk surface, here at relative heights z/r = 0.4-0.5. In this re- 
gion, PAH heating is usually the most important heating process. 
Therefore, this line shows a strong correlation with the assumed 
PAH abundance, /pah, UV irradiation, and disk flaring. 

The optical [OI]6300A line originates from hot gas in the 
inner disk, about 0.06- 0.7 AU in this model. As a forbidden 
line, with very small Einstein coefficient A41 ~ 6.5 x 10~ 3 s _1 , 



5 For more typically extended disks, R out = 100 - 500 AU, the 
[OI] 63 [im line mo stly originates in 10 - 100 AU, slightly larger in case 
of Herbig Ae disks dKamp et alll2010h . 
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it is mostly optically thin Tii ne ~ 0.1-1. With an excitation 
energy of ~ 15900 K, it probes the existence of hot and dense 
gas T gas ^ 3000 K in front of an optically thick continuum. To 
our surprise, this line is not significantly excited by OH photo- 
dissociation in the model. Neglecting this pumping effect (see 
Appendix lA.5b weakens the [OI] 6300 A line flux by only 14%. 
Thus, the [OI] 6300 A line comes from the bottom of the hot 
atomic layer, from about z/r = 0.3-0.4 in this model. As soon as 
molecules form, for example OH and CO, the temperature drops 
significantly by molecular line cooling, and the [OI] 6300 A line 
cannot be excited any longer. Since the line is optically thin, its 
flux reacts quite sensitively on disk mass. We note that the pre- 
dicted line FWHM = 35.7 km/s is in excellent agreement with 
the observations, FWHM = (38+ 15) km/s, meaning that the spa- 
tial origin of the [OI] 6300 A line is about correct in the model. 

The near-IR o-H 2 v=1^0S(1)2.122fim line also probes 
the inner disk regions, 0.05- 0.6 AU in this model. The line 
is quite optically thick, Ti; ne as 5-10 in this region, where the 
dust is also optically thick. Since its Einstein coefficient is ex- 
tremely small, A$i *7x 10~ 7 s _1 , the line needs large H2 column 
densities of order 10 23 cm~ 2 to become visible over the con- 
tinuum level. Such column densities are only reached in quite 
deep layers, z/r w 0.15, so this line forms in much deeper lay- 
ers than the [OI] 6300 A line. The temperature contrast between 
gas and dust is already quite small in these layers (Fig. IC.21 . 
which limits the 0-H2 2. 122 /mi line flux. The excitation energy 
of the line is about 7000 K, i.e. r gas <: 1000 K is required to 
collisionally excite it. In the best-fitting model such conditions, 
Tgas ^ 1000K» 7dust, are provided by the effect of exothermic 
chemical reactions (see Appendix IA.8I 1 which is active in the 
warm and dense gas close to the inner rim, deep but not too deep, 
just where this line forms. Therefore, we see a clear correlation 
between the 0-H2 2. 122 yizm line flux and the assumed heating ef- 
ficiency of exothermic reactions y chem . The line is also substan- 
tially pumped by ^-formation on grains (see Appendix IA.61 . 
When neglecting this effect, the line attains a flux of only 35% 
of the value from the full model. We conclude that the H2 line 
is sensitive to a temperature contrast between gas and dust, and 
H2-formation, in quite deep disk layers close to the inner rim. 
However, the predicted line is too broad (FWHM = 37 km/s) as 
compared to the observations (FWHM = 18 + 1.2 km/s). We were 
unable to find any disk model, among the ~ 20000 models com- 
puted, that shows such a narrow o-H2 2.122jt/m line. Thus, the 
line forms too close to the star in the model, and we must be 
very careful when drawing our conclusions about the nature of 
ET Cha from the observed 0-H2 2. 122 //m line. 

The sub-mm CO J = 3 -> 2 line (not depicted in Fig. Qj} 
probes the size of the disk and the gas temperature in the out- 
ermost disk layers. It is massively optically thick, rii ne <: 10 5 in 
the line forming region 2-7 AU, where the dust is optically thin. 
Its flux is roughly proportional to the projected disk area times 
the gas temperature at relative heights z/r x 0.4-0.5 in these 
outermost parts of the disk. 

5.6. Disk gas mass and disk shape 

The determination of the total disk gas mass of ETCha in 
this paper is based on three detected gas emission lines, 
namely [OI]63//m, [OI]6300A (LVC), and o-H 2 v=1^0 
S(l)2.122/im, which probe complementary radial and vertical 
disk regions. However, the first line, [OI] 63 /im, is massively 
optically thick, and the latter two lines mainly probe the hot gas 
in the inner disk regions. CO J = 3 — > 2 was not detected. These 



circumstances immediately suggest that a robust gas mass deter- 
mination of ET Cha is difficult. 

Table [7] shows that various solutions in form of well-fitting 
disk models can be found where the total gas mass differs by up 
to two orders of magnitudes (factor 300). A careless reading of 
Table|7]would suggest that the particular values for M gas , e and /3 
do not matter much, but this is not true. Careful inspection of the 
solutions in Table |7]by systematic variation of selected parame- 
ters (see e. g. Fig. [9] for the best-fitting model) shows that these 
solutions all represent well-defined local minima, where small 
changes of any parameter leads to a considerable deterioration 
of the combined line + continuum fit quality x- If one would 
measure the uncertainty of gas mass determination from these 
dependencies alone (for instance the deviation AM gas where % 
doubles) one would arrive at relatively small errors, about 40%. 
But such an error estimate would be misleading as it would not 
take into account the complicated manifold of local minima in 
parameter space. 

All solutions in Table [7] fit our entire set of line and contin- 
uum observations about equally well. The different values for 
the gas mass, however, come in certain combinations with other 
disk shape parameters, like the column density power-law expo- 
nent e and the flaring power f3. Certain fine-tuned combinations 
of M gas , £ and /3 do apparently all provide the proper density and 
temperature conditions in the disk that result in almost exactly 
the same observables, for example small M gas in combination 
with large e and large /3, or vice versa. 

These relations can be understood by minimum column den- 
sities of warm H2 and atomic oxygen in the inner disk parts 
that are inevitably required to make the [OI] 6300 A (LVC) and 
o-H2 2.122/im lines visible over the strong continuum at opti- 
cal and near IR wavelengths. For example, in the least massive 
model 1 in Table [7] there is so little gas in the disk that a rela- 
tively large value of e- 1.16 is required to concentrate the mass 
in the inner disk parts and so to provide the necessary column 
densities there. All line fluxes are in good agreement with the 
observations then, but the model has problems to find a good fit 
of the 10/im and 20 pm silicate features. Also the flaring index 
of p = 1.33 is quite extreme, leading to a relative disk height of 
about z/r = 1 .2 in the outer disk parts - such a "disk" would be 
taller than wide. From these arguments, we will discard model 1 
from our selection of valid solutions in the following, and claim 
a minimum gas mass of about 5 X 10 -5 M Q to achieve the nec- 
essary minimum gas column densities in the inner disk regions 
and to fit the silicate dust emission features simultaneously. 

Concerning the other direction, we find it hard to provide any 
solid argument why the disk gas mass must not exceed a cer- 
tain maximum value. In combination with our quite robust dust 
mass determination of Md us t ~ 3 x 10~ 8 M Q , however, even our 
minimum gas mass of 5 x 10~ 5 M o already implies a very high 
gas/dust ratio of 2000. An even higher gas mass would translate 
into an accordingly higher gas/dust ratio. The largest value we 
found with our evolutionary strategy is 2.5 x 10~ 3 M . 

Summarising our modelling efforts, after having tried about 
20000 disk models and various choices of input physics, we con- 
clude that the disk gas mass of ET Cha is only little constrained, 
our confidence interval is about M gas » (5 x 10~ 5 - 3 x 10~ 3 ) M Q . 
The surface density exponent e and the flaring parameter /3 are 
likewise poorly constrained, although the models favour small 
e^0 to fit the silicate dust emission features. 

In view of the outflow discussion (Sect. 16.41 ). we note that 
the unclear physical origin of [OI] 63 fim obviously increases 
our uncertainty in gas mass determination. We have set up a 
similar evolutionary optimisation run as depicted in Fig. [8] but 
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now treating the [OI] 63 /urn line flux, as emitted by the disk, 
as an upper limit with < 2 x 1CT 18 W/m 2 . This run resulted in 
a disk of similar mass as compared to the best-fitting model, 
M gas = 5.5xl(T 4 M , but of even smaller size R out = 2.1 AU, with 
an [OI] 63 /mi line flux of ~ 4 x 1(T 18 W/m 2 . We conclude that 
the [OI] 63 /im emission line flux is not crucial for the disk gas 
mass determination of ET Cha, because this disk is tiny and the 
hot gas is responsible for our detected [OI] 6300 A (LVC) and 
o-H2 2.122yum lines is located in the inner regions. For other, 
more extended protoplanetary disks, the [OI] 63 /.im line will be 
less optically thick and hence more useful for the purpose of 
gas mass de t ermin ation dKamp et al.L 120101 : iPinte et all [2010; 
IWoitkeetal.Ll20lol) . 

6. Discussion 

6.1. Scale heights and unidentified heating 

From the SED fitting, we get surprisingly robust results con- 
cerning the assumed vertical disk scale height, Hq = (0.008 - 
0.01 1) AU at reference radius r = 0.1 AU, throughout all SED- 
fitting models. In the model, the reference scale height in combi- 
nation with the flaring power (3 determines how much star light 
is captured by the disk. This light (unless scattered away) is ab- 
sorbed by the dust in the disk surface and then thermally re- 
emitted, heating also the inner disk parts. Thus, Hq regulates the 
dust temperature in the disk. Since ET Cha is very bright in the 
(3 - 8) pm region, we require large relative scale heights of the 
order of 10% close to the star to create enough dust emission 
from warm and hot grains to fit the near and mid IR. 

Appendix lD.4l (see lower left plot in Fig. IC.41 > demonstrates, 
however, that our prescribed scale heights are significantly 
larger, by a factor of 2-3, than those derived from self-consistent 
models where vertical hydrostatic equilibrium is assumed. This 
mismatch can be interpreted in three ways: (1) the close mid- 
plane regions of T Tauri disks are not in hydrostatic equilibrium, 
(2) the midplane temperatures are actually 4-9 times higher than 
assumed in the model (assuming H oc cj oc yT), or (3) there is 
an additional dust heating process active in the close disk mid- 
plane regions that leads to an additional energy flux through the 
dust component resulting in more observable near-mid IR pho- 
tons without changing the temperatures much. Possibility 1 can- 
not be excluded per se. Possibility 2 seems unrealistic because 
it would cause the dust to evaporate. We favour possibility 3. 
Appendix ID. 51 shows that the inclusion of non-radiative dust 
heating via inelastic gas-dust collisions (thermal accommoda- 
tion, driven by gas-dust temperature differences created through 
exothermic chemical reactions) leads to similar effects than in- 
creasing the scale height. Because of the p 2 -scaling of chemical 
reactions, this additional heating affects the dust preferentially in 
high-density regions, increasing the production of near IR pho- 
tons needed to fit the SED of ET Cha with smaller scale heights. 
By means of an extra run of the evolutionary strategy, we found 
out that we can reduce the scale height by about 35% to fit the 
SED, if we include this effect. 

More unidentified dust heating processes may be active in 
the close midplanes of T Tauri dis ks. However, v iscous heating 
(according to the formulation bv lFrank et all [T992h is not doing 
a particularly good job in explaining the scale height inconsis- 
tencies of ET Cha, see Appendix lD.3l Viscous heating dumps ad- 
ditional energy ocp, i.e., after volume-integration, preferentially 
into the outer cold regions, leading to the production of more far 
IR photons. Moreover, viscous heating simply has little effect on 
the resulting disk temperatures and observable continuum flux in 



case of low-mass disks as for ET C ha. In fact, the vis cous heat- 
ing according to the formulation by dFrank et al.l [T992) produces 
artifacts in the uppermost tenuous disk layers, where the viscous 
heating ocp cannot be balanced by any cooling ocp 2 . 

6.2. Effects of other physical processes 

Appendix [D] discusses a number of further physical processes 
and their influence on the model results that have not been men- 
tioned so far and have not been included in the models discussed 
so far. To summarise, we find that 

(i) X-rays have very little effect (see Appendix ID. 2\ and the 
neglecting of X-rays in our main model for ET Cha is fully jus- 
tified. An X -ray luminosity of Lx = 6 x 10 28 erg/s, as observed 
for ETCha (lLopez- Santiago et all 1201 Ol) . turns out to be much 
less important for the disk heating and ionisation as compared to 
the strong FUV of ETCha, Luv = 6.5 x 10 31 erg/s as measured 
between 912 and 2500A, based on our Hst/Cos and Hst/Stis 
observations. 

(ii) The treatment of H2-formation on grain surfaces is one of 
the most important yet quite uncertain processes for astrochemi- 
cal modelling. Appendix lD.6l shows that different approaches to 
calculate the ^-formation rate lead to a systematic uncertainty 
in the model for the computation of the o-H2 2.122yum about a 
factor of 5. Other spectral lines are less effected. 

6.3. Disk inclination and outflow velocity 

The low blue-shifts (-42 km/s) seen in the optical emission lines 
(Fig.O are quite unusual for the jets/outfl ows of T Tauri stars . 
Typical outflow velocities are ^ 100 km/s (Harti gaii et al.L fl995). 
A possible explanation could be the projection effect in case of a 
strongly inclined disk (i.e., close to edge-on). Adopting 100 km/s 
as a lower limit for the outflow velocity of ET Cha would suggest 
a disk inclination of >65°. 

However, as argued in Sect. 14.31 we find that inclinations 
in the range 0°-50° are in agreement with the observations, but 
larger inclinations would result in a partial obscuration of the star 
by the disk, with dramatic effects on the SED. The SED analy- 
sis therefore leads to the conclusion that the inclination of the 
system, as measured from face-on, is 50° or less which, in turn, 
translates to an outflow velocity of 65 km/s at most. If this inter- 
pretation is correct, this would make the outflow from ET Cha 
one of the slowest known outflows from a T Tauri star. 

6.4. Outflow and disk lifetime 

Our analysis of the optical [OI] 6300 A and [SII] 6731 A emis- 
sion line profiles and fluxes (see Appendix 151 results in an es- 
timate of the outflow mass-loss rate of ETCha of M out g ow * 
10~ 9 M o /yr. Such an outflow can contribute to the [OI]63pm 
emission, which would render the [OI] 63 //m line flux, as emit- 
ted from the disk, smaller as assumed in Sect. 14.51 In summary, 
Appendix [B] shows that both app roaches, the simple energetic 
outflow analysis by iHollenbachl (Il985[). a s well as the shock 
models computed by dHartigan et alll2004l) . suggest that the out- 
flow from ETCha does contribute a substantial, if not domi- 
nant, fraction of the observed [OI] 63.2pm emission line flux. 
However, without spatially or velocity resolved observations of 
the outflow, and its subsequent detailed modelling (rather than 
using some "template" shock models), it is presently impossible 
to assess the exact contribution of the outflow to the line emis- 
sion. 
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Our estimate of the outflow rate of ET Cha is large compared 
to the total disk mass we derive, M^k ~ 6 x 1CT 4 M , as it sug- 
gests a disk lifetime of only M gas /M out fl OW * 0.7 Myr, inconsis- 
tent with the age of the Chamaeleontis cluster of ~ 8 Myr. If we 
would assume a generic gas/dust ratio of 100 and take our dust 
mass determination for granted, Md ust as 3 x 10~ 8 M , it is even 
worse. The disk lifetime would be even shorter in this case, only 
3000 yrs, way too short to be feasible, unless we are just observ- 
ing a temporary but short-lived lOOx peak in outflow rate. 

We also notice that the mass accretion rate of ET Cha was 
estimated by lLawson et al.l d2004t) to be equally large, M acc as 
10~ 9 M o /yr as the outflow mass loss rate, leading to simi- 
lar lifetime inconsistencies. Furthermore, a branching ratio of 
^outflow / M acc ~ 1 is highly unusual, a few perce nt seems to be a 
well-e stablished value for T Tauri stars (see e.g. Hartiga n et al.L 
119951). 

iMurphv et al] (1201 ll) reported on highly variable H a equiv- 
alent widths and, accordingly, mass accretion rate, for the old 
T Tauri stars in the 77 Cha cluster. A factor of 100 variation in 
the accretion rate is observed in one newly-identified halo mem- 
ber of 77 Cha. But only a simultaneous reduction of both M acc 
and M out fl ow would help to resolve the aforementioned lifetime 
inconsistency, which doesn't seem very likely. 

One possible explanation would be a massive but short-lived 
outflow due to a flare from a former epoch, when the mass ac- 
cretion rate was at least 10 times higher. An outflow of lOOkm/s 
would need ~ 50 yrs to reach a distance of 1000 AU. This would 
be in agreement with the somewhat slow outflow velocity we de- 
rive, 42-60km/s, because the outflow might have slowed down 
ever since. However, the optical emission lines do not show any 
evidence for a red-shifted HVC, as one would expect for a sym- 
metric bi-polar outflow. The only plausible explanation of the 
missing red-shifted HVCs is that these components from the far 
side are attenuated by the dust in the disk. According to our disk 
model, the disk is optically thick at 6300 A up to the outer radius 
~ 10 AU (see Fig.[TTJ. Therefore, the line emitting region of the 
blue-shifted optical emission lines must be quite small, less than 
5 AU if seen under 60° disk inclination, which corresponds to an 
age of only 0.25 yr. 

7. Summary and conclusions 

This paper has reported on new observations of ET Cha with sev- 
eral instruments: Herschel/Pasc, Ctio/Andicam, Hst/Cos/Stis, 
and Apex. In combination with published data from Spitzer, 
Gemini/Phoenix and Aat/Ucles, we have collected an unprece- 
dented observational data set about this object, including pho- 
tometry, UV spectra, high-resolution optical spectrum, near and 
mid IR spectra, and far IR and sub-mm line fluxes. 

We have calculated united gas and dust models for the disk 
of ETCha that can simultaneously fit all line and continuum ob- 
servations, except for a too broad o-H2 2.122yum emission line 
profile. The observations also show some blue-shifted compo- 
nents of optical emission lines that point to an outflow and are 
not included in the models. 

This paper has explored the parameter space of the disk mod- 
els by using an evolutionary strategy to minimise the discrepan- 
cies between model predictions and observations. The paper has 
also introduced a number of basic improvements to the ProDiMo 
disk modelling code concerning the treatment of PAH ionisation 
balance and heating, heating by exothermic chemical reactions, 
several non-thermal pumping mechanisms for selected gas emis- 
sion lines, and formal solutions of the line transfer problem at 
given inclination (AppendixlAl. 



From the disk modelling we find a rich variety of fitting disk 
models that can explain our observations about equally well. 
Some of the model parameters (like the dust mass and the outer 
radius) can be determined with some confidence whereas other 
parameters (like the disk gas mass) are poorly constrained: 

- The new Herschel/Pacs photometric fluxes at 70 fim and 
160^/m constrain the disk dust mass of ETCha to be about 
.Mdust = (2-5) X 10~ 8 M Q , putting the object at the borderline 
between optically thin and optically thick. 

- Then strong near IR excess of ETCha can be fitted 
with a disk that is truncated at R[ n as 0.022 AU (where 
Tdust ~ 1500 K) which is located slightly outside of the 
co-rotation radius of 0.015 AU. The latter is calculated 
according to the assumption that the star rotates with a 
period P = 1.7 days as suggested by our v lot sin(z') analysis 
of rotationally broadened stellar absorption lines. 

- From the Apex CO J — 3 — > 2 non-detection, we can infer, 
with confidence, that the disk of ETCha must be tiny in 
radius. The models favour an outer disk radius as small 
as R oia as (6 - 9) AU. All disk models with R om ^ 25 AU 
would violate the 3cr CO J — 3 — > 2 non-detection limit, 
independent of chemical details. 

- The SED-fitting suggests that the dust grains in the surface 
of the protoplanetary disk of ETCha (where the near-mid 
IR continuum forms) must be small in radius (sub-micron 
sized) and opaque in the optical and near- mid IR, i. e. 
absorption must dominate over scattering opacities. In 
the models, these spectral properties are provided by the 
inclusion of about 25% amorphous carbon, but other options 
like metallic iron are also possible. 

- The disk gas mass of ET Cha is poorly constrained by our 
line observations CO J — 3 — > 2 (non-detection), [OI] 63 yum, 
[OI]6300A (LVC), and o-H 2 v=1^0S(l)2.122yum. We 
find a variety of about equally well fitting disk models with 
total gas masses M gas = (5 x 10~ 5 - 3 x 1O~ 3 )M . The 
forbidden lines of [OI]6300A (LVC), and o-H 2 2.122//m 
are close to optically thin and hence quite useful for gas 
mass determination, but originate from hot gas in the inner 
disk regions only. [OI] 63 /im is massively optically thick 
and does hence not discriminate much between the various 
disk models for the case of ETCha. We would need to 
observe additional far-IR or sub-mm spectral lines that 
originate from the outer disk parts to determine the gas 
mass with more confidence. However, in order to explain 
the [OI] 6300 A (LVC) and o-H 2 2.122 /an lines with a disk 
model, we require line optical depths Ti; ne ^ 1 which trans- 
lates into total vertical column densities in the inner disk 
regions of about No ~ 10 21 crrT 2 and Nm ^ 10 24 crrT 2 for 
atomic oxygen and molecular hydrogen, respectively. Such 
column densities are incompatible with M gas <5 x 10~ 5 M Q . 

- The wide range of fitting gas masses is related to particular 
values of the disk shape parameters, namely the column 
density exponent e and the disk flaring power /3. Large gas 
masses need to be combined with small values for e and /3, 
and vice versa. 

- From our SED modelling of ETCha we derive disk scale 
heights of the order of 10% relative to radius close to the star, 
which is about a factor of 2-3 larger than the scale heights 
inferred form self-consistent models that assume hydrostatic 
equilibrium. This discrepancy can partly be explained by an 
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additional non-radiative heating of the dust close to the star, 
for example via exothermic chemical reactions. 

These results suggest a surprisingly high value for the overall 
gas/dust mass ratio of ET Cha of at least 2000, or even 20000. 
Whether or not the gas responsible for the [OI] 6300 A and o- 
H2 2.122/im emissions is still physically connected to the disk 
is not entirely clear, but the observations show that this gas is 
at least not moving much with respect to the disk. Possibly, the 
overwhelming majority of dust particles responsible for the near 
to far-IR emission of the star has already been transformed into 
larger pebbles or solid bodies, which have negligible opacities. 

The fluxes in blue-shifted emission line components like 
[OI]6300A (HVC) suggest an outflow with mass-loss rate 
^outflow ~ 10~ 9 M /yr, on a similar level as the reported mass ac- 
cretion rate of ET Cha. The low velocities (=s -42 km/s) suggests 
a high inclination angle of at least 60° (rather more). However, 
inclinations in excess of 50° are inconsistent with our SED- 
modelling, which favours i ^ 40°, from which we determine an 
outflow velocity of 65 km/s at most (typical values ar e in excess 
of 1 00 km/s for T Tauri outflows lHartigan et all 1 19951) . If this in- 
terpretation is correct, ET Cha would possess one of the slowest 
known outflows from a T Tauri star. 

An outflow with mass loss rate M ollt fl ow * 10~ 9 M o /yr is 
likely to contribute significantly to the [OI] 63/mi line flux as 
observed with Herschel/Pacs. Given the observational data we 
have collected in this paper, we cannot discriminate between out- 
flow or disk origin of [OI] 63/im. 

To conclude, ET Cha seems to be an extraordinary and puz- 
zling object concerning the evolution of protoplanetary disks. 
Despite its age of about (6 - 8) Myr, there is evidence of active 
accretion onto the central star, on a similar level as the derived 
outflow mass-loss rate. According to the low disk masses we de- 
rive, the disk lifetime M gas /M acc as (0.05 -3) Myr is inconsistent 
with cluster age. If a generic gas/dust ratio of 100 was assumed, 
the disk lifetime, based on our relatively robust determination 
of the disk dust mass, would be even shorter, only ~ 3000 yrs. 
Either the object is actually much younger than the age of the 
Chamaeleontis cluster or the object is going through a phase of 
unusually high mass accretion rate and outflow. 
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Appendix A: New features in ProDiMo 

A. 1 . Fixed disk structure 

In contrast to earlier publications, we use a parametrised descrip- 
tion for the shape and distribution of gas and dust in the disk in 
this paper 



p(r, z) = po 



(-) 

\r f 



: H 
H(r) 



exp 



( 2H(r) 2 ) 



(A.l) 



between an inner and outer disk radius, R m and R out , respectively, 
with sharp edges. p(r, z) is the local gas mass density. H(r) is the 
vertical scale height of the disk, assuming to vary with radius as 



H(r) = H (-) 
\r / 



(A.2) 



Hq is the reference scale height at reference radius ro- e is 
the column density power-law index and ft the flaring power. 
The constant po is adjusted such that the integrated disk mass 
2k ffp(x,z)dz rdr equals Mdisk- 

A.2. Dust size distribution and dust settling 

The dust grains are assumed to have a power-law size distribu- 
tion in the unsettled case as 



/(a) oc a ' 



(A3) 



between minimum and maximum grain radius, a m i n and a max , 
respectively. The free constant in Eq. (lA.3t is adjusted to result 
in the prescribed unsettled dust/gas mass ratio pa/p. 

A very simple recipe has been implemented to account for 
the major effects of vertical dust settling. We assume that the 
dust grains are distributed vertically with a smaller scale height 



H'(a, r) — H(r) ■ max Jl, a/a s 



-s/2 



(A.4) 



where H(r) is the gas scale height, and s and a s are two free 
parameters. Since ProDiMo can self-consistently calculate the 
vertical stratification of gas from the resulting gas tempera- 
tures and mean molecular weights, in which case H(r) does not 
exist, Eq. dA.4l > can not be used directly in the general case. 
Instead, we make use of the equation that defines H(r), namely 
CT=H(r) £2(r) and write 



c' T (a) - cj ■ max 



{l,a/fl s 



-s/2 



(A.5) 



where cj is the local gas isothermal sound speed, c' T (a) is the re- 
duced variant for dust size a, Q. = VKepier/ r is the angular velocity 
and VKepier the Keplerian velocity. Equation JA.Sb is then used 
to calculate the vertical distribution of dust particles of differ- 
ent sizes with respect to the already determined gas stratification 
by treating them like an independent fluids with smaller cj as 
compared to the gas 



1 dp'(a) 
p'(a) dz 



Idp 

p dz 



(A.6) 



where p — c\ p belongs to the gas, and p'(a)-Cj(a) p'(a) belongs 
to dust particles in [a, a + da]. g z is the local z-component of 
gravity. The solution of the differential Eq. ( IA.61 > is 



f(a,z) 



const(a) 

c'j(a,z) 



exp 



'|^ ta (4«/ ( »,z,) 



(A.7) 



where f'(a,z)da = p' (a) I m&{a) is the settled dust size distri- 
bution function, and ct(z) and f(a,z) oc p(z) are considered as 
known. The constant const(a) is then determined to make sure 
that the vertical column density of dust particles for every size is 
conserved 



jf'(a,z)dz= Jf(a,z)dz. 



(A.8) 



According to these assumptions, all dust quantities that used to 
be constant, such as the local dust/gas mass ratio pd/p, the dust 
moments (a), (a 2 ), (a 3 ), and the dust opacities per mass, be- 
come spatially dependent quantities. The dust absorption and 
scattering opacities are calculated by applying effective mixing 
dBruggemanlll935l) and Mie theory, based on f'(a, z). 

A.3. PAH ionisation equilibrium and PAH-heating: 

We consider a typical size of PAH molecules with Nq = 
54 carbon atoms and A^h = 18 hydrogen atoms (circum- 
coronene), motivated by studies that PAHs with much less car- 
bon atoms would not be stable around young stars on timescales 
of a few Myr. Much larger PAHs are not consistent with 
the spati al extent of o bserve d PAH emission in various bands 
(see e.g. IVisser et all 120071) . We include PAH , PAH, PAH + , 
PAH 2+ and PAH 3+ as additional specimen in the chemical re- 
action network. Circumcoronene is probably among the small- 
est PAHs th a t can survive in a disk around Herbig Ae stars 
(IVisser et all 120071) . The follo wing processe s are considered 
in detail: PAH-photoionisation dTielensl 120051) . el ectron recom- 
binati on and some charge exchange rea ctions dWolfire et all 
120081 E ower & Pineau des Foret si l200l . These reactions do 
not change the basic PAH lattice, but only affect the charging 
of the PAH molecules. Hence, the total amount of PAHs is con- 
served and treated like an element with given element abundance 



e(PAH)=/ PAH <™ 



(A.9) 



where XfJ^ = 3 x 1(T 7 being the sta ndard ISM par ticle abun- 
dance with respect to hydrogen nuclei (Tielens, 2008), and /pah 
is the fraction thereof assumed to be present in the disk. 

Concerning the photo-ionisation reactions, we consider the 
PAH absorption cross sections o"p AH (v) [cm 2 ] from lLi & Drainel 



(2001, see their Eqs. 6-12 and A2-A3) with rec ent updates for 
the resonance parameters from lDraine &Lil(l2007l) . We use their 
neutral PAH cross section for charge k = 0, and the charged 
cross section otherwise. The photo-ionisation rates [s _1 ] for PAH 
molecules with charge k are calculated as 



- f Mhr 

he J912A 



°pah<» v/ v (r,z) Y k v s v (r,z) dA 



(A. 10) 



where v[Hz] is the frequency, J v (r, z) [erg/cm 2 /s/Hz/sr] the 
mean intensity as computed by the dust continuum radiative 
transfer, and s v (r, z) the following self-shielding factor 



r P v AH (r,z) 



e(PAH) < ( r)p AH (v) Min [N™(r, z), N^(r, z 



s v (r,z) = exp(- T P AH (r,z)) 



)}(A.ll) 
(A. 12) 



and N r ^(r,z) are the radial inward and vertical upward 

hydrogen nuclei column densities [cm -2 ] in the disk, as mea- 
sured from point (r, z). For simplicity, we put (cr) PAH (y) = 
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Table A.l . Larger non-LTE m odel atom s. Reference s are a s fol- 



lows: A =ISchoier et al l (120051). NIST =lRalchenkol d2009h K 



iKrems et alJd2006l) . S = IStorzer & H ollenbachl d20o' 



species 


#levels #lines 


coll. partners 


references 


OI 


91 647 


p-H 2 ,o-H 2 ,H,H + ,e- 


A, NIST, K, S 


CI 


59 117 


p-H 2 ,o-H 2 ,H,H + ,He,e- 


A, NIST 


CII 


10 31 


p-H 2 ,o-H 2 , lie- 


A, NIST 



0.5 (o"p AH (v) + (Tp^ H ( v)). The photo-electron yield is taken from 
dJochimset all 1 19961) 



the three spin-orbit split ground states 3 P exist with electrons and 
hydrogen atoms as coll ision partners dStorzer & Hollen bach. 
2000; IKrems et all 2006). For C and C + , we use the ra t es col - 
lected in the Leiden Lambda database (Schoi er" et all [2005), 
whose collision rates are mostly limited to the ground state elec- 
tronic levels. 



A.5. [OI]6300A -pumping by OH photo-dissociation 



Y k = 



1 

hv-IP 



9.2 eV 




, hv > IP k + 9.2eV 

, IP k + 9.2eV >hv> IP k 

, hv < IP k 



Photo-dissociation of OH by absorption of UV photons into the 
1 2 lT state, A ^ (973—1 350) A, produces an electronically excited 
oxygen atom O('Z)) as 



(A.13) OH + hv— > 0('D) + H, 



(A. 16) 



According to Weingartn er & Drainel (1200 lb . the ionisation po- 
tentials IP k depend on PAH-size and charge k. Using their 
Eqs.(l) and (2) for N c = 54, the results are IP 1 = 3.10eV, 
IP = 6.24eV, IP +l = 9.38eV and IP +2 = 12.5eV. The for- 
mula is an approximation that reproduces the ionisation potential 
of benzene as the smallest PAH and graphite as infinitely large 
PAH. The threshold wavelengths are given by AjL = hc/IP k . 

With these photo-ionisation and recombination rates, the lo- 
cal particle densities of PAH , PAH, PAH + , PAH 2+ and PAH 3+ 
are calculated consistently with the gas-phase and ice chemistry, Cf^ m = 0.55 x «oh / «o 
with a considerable influence on the resultant local electron den- 



while absorption into the 1 'A and 3 2 Y1 states, A ^ 
(1100- 1907) A, leads to the g round state 3 P. According to 
van Dishoeck & Dalgarnol dl984l) . about 55% of the OH pho- 
todissociations by an interstellar UV fielc@ result in an O-atom 
in the l D 2 state, which happens to be the upper level of the 
63 00 A line l D 2 — > ^P 2 - This chemical pumping was suggested 
bv lAcke et al.1 d2005l) to play a major role in the formation of the 
6300 A line around Herbig Ae/Be stars. We take this effect into 
account by introducing a quasi-collisional pumping rate [s~'] 



(A. 17) 



sity. Once these particle densities «p AH have been determined, 
the total PAH heating rate [erg/cm 3 /s] can be calculated as 



PAH 



An 

he 



"ihi 

X"pah f o- k mn (v)vJ v Y k s v (hv-IP k )dA (A. 14) 

I, ^ 



912A 



The PAH recombination cooling rate is calculated as 



A 



' 1 PAH ' 



' ^PAH^g) ( 1 -5 Wg) 



where R^ 1 is the OH photodissociation rate [1/s]. We apply 

Eq. $AT?} to / = {1,2,3}, where 1 = 3 P 2 , 2 = 3 Pi, 3 = 3 P , and« = 
4 = l D 2 . We add these excitation ra tes to the other col lisional ex- 
citation rates C; M (see Sect. 6.1 in Woit ke et all 12009). In writing 
Eq. dA.17l >. we rely on statistical equilibrium and assume that the 
majority of chemical reaction channels forming OH start from 
the 3 lowest 3 P levels with no ~no(l)+"o(2)+«o(3), making sure 

tVmt n I "oW /^chem , «o(2) / ^ c hem , «o(3) ^cheiiA _ n cc nOH 

that n y— C I4 +— C 24 +— C 34 j - 0.55 «OH^ ph , 

no© 



(A. 15) where -2" are the fractional populations. 



where ^ PAH (r g ) 



cm 3 s 1 ] is the PAH recombination rate coeffi- 



cient, n e is the electron density and the gas temperature. As a 
general result, we obtain a top-down layered PAH charge struc- 
ture, with PAH 2+ in the uppermost layers with reduced PAH- 
heating, but PAH" and intensified PAH-heating close to the mid- 
plane. For massive disks (unlike ETCha), where the electron 
density virtually vanishes in the deepest layers, the PAH charge 
becomes neutral again in the midplane. 

A.4. Fluorescent UV-pumping 

We have replaced th e small model atoms for O, C and C + (as 
listed in Table 4 of Woitk e et all 120091) by larger ones, see 
Table I A. 1 1 (this paper), to account for fluorescent UV and op- 
tical pumping of the lower levels responsible for the far IR fine- 
structure lines. The new model atoms have collisional data only 
among the lowest states, but much more radiative data. The ad- 
ditional energy levels and transition probabilities are taken from 
the National Institu te of Standards an d Technology atomic spec- 
troscopic database ( Ralc henkol l2009t) . We selected all evaluated 
transitions longwardof 912 A. Rovibronic transitions arepermit- 
ted and their probabilities are of the order of (1 - 10 6 ) s _1 , much 
higher than typical collision rates (~10~ 9 s -1 ). 

Collisional data connecting the two first electronic-excited 
states l D 2 at 22830 K and l S at 48370 K of neutral oxygen and 



A.6. H 2 -pumping by its formation on dust surfaces 

The formation of H2 on grain surfaces liberates the binding en- 
ergy, causing the newly formed H2 mol ecules to leave the s urface 
in a vibrationally highly excited state. Dulev & Williams (1986) 
estimated that the H2 molecules will be released in vibrational 
states as high as v$7, but in the rotational ground state J-0 for 
P-H2 and 7=1 for 0-H2. 



H + H + dust H 2 (v<7,J = 0/l) + dust 
We treat this formation-pumping by 

C^ m =«*"H^ t /"to, 



(A.18) 



(A. 19) 



where R d ^ t is the H2 formation rate on dust surfaces [1/s], u is 
the index of the highest state with v < 7 and J < 1 that still has 
collisional de-excitation rates in the database, and Eq. ( IA.19I ) 
is applied to all states I < u, i. e. to all states that have lower 
excitation energy than u. n lot is the total ortho or para H2 par- 
ticle density, respectively, and a = n /p-n 2 /nn 2 is the prescribed 
ortho/para-H2 fraction. 



6 We note that the branching ratio could be different for our adopted 
star+UV spectrum. A maximum branching ratio of 100% would would 
increase our [OI] 6300 A line predictions by less than a factor of 1.8. 
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Our formulation of the H2-pumping, therefore, leads to a 
constant de-population of all states I < u with a certain timescale 
C* em [s -1 ] (simulating the disappearance of H2 in all states due 
to other chemical channels that are feeding the H2 formation) 
and an equally large population [cm 4 s _1 ] of the state u due to 
the formation of H2 on grains. 

The condition of state u being collisionally connected avoids 
artifacts at very low densities, where collisions are rare and 
the pumping would lead to an almost complete de-population 
of all low-lying states. According to our current collection of 
collisional H2 de-excitation rates, this condition implies u - 
(v = 3,7 = 0/1). We checked that our modelling of the o- 
H2V=1— >0S(1) line at 2.122jum does not depend much on the 
choice of this upper level u as long as v > 1 . 



A. 7. Line transfer 

During solving the chemistry and energy balance of the gas, 
ProDiMo calculates the various level populations of atoms, ions 
and molecules by means o f an es cape probability method (see 
Sect. 6.1.1 in [Woitke et al.L 120091) . These are stored and used 
later to perform a formal solution of line transfer for selected 
spectral lines. 

The continuum + line radiative transfer equation is given by 



dl 



(A.20) 



where l v is the spectral intensity and s the distance on a ray. The 
source function and extinction coefficient are given by 



= £ 



+ 0^ • 



(A.21) 
(A.22) 



Assuming isotropic dust scattering, the continuum and line 
transfer coefficients are given by 



D _ abs sea 



t _ hv 
£^ — ~. n uA u l 
An 

hv 



1 " v I \ 
<T = —\niB h , - n u B u ij , 



(A.23) 
(A.24) 

(A.25) 
(A.26) 



where /cj bs and ^ ca [cm -1 ] are the local dust absorption and scat- 
tering coefficients, B Y is the Planck function, J v is the local mean 
intensity (taken from the results of the continuum radiative trans- 
fer), v is the line centre frequency, n u and n\ are the level popu- 
lations [cm 4 ] in the upper and lower state, respectively, and A u \ 
and B„i, Bi u are the Einstein coefficients. 

The profile function <p v [Hz -1 ] is assumed to be given by a 
Gaussian with thermal + turbulent broadening 



v + n ■ v 

c 



v y/nAv ^ \ Av 2 



(A.27) 
(A.28) 



where v' is the observers velocity with respect to the star along 
backward ray direction n, v is the 3D-velocity of the emitting 
gas with respect to the star (assumed to be given by Keplerian 
orbits), and x the local velocity shift. The velocity width is Av 2 = 
vl + v 2 urb and the thermal velocity v 2 h = 2kT g /m where m is the 
mass of the line emitting species. 



We u s e the same setup of parallel rays as described in 
iThi et alJ (l2010bl see their Sect. 2.3), organised in about 150 
log-equidistant concentric rings in the image plane, each sub- 
divided into 72 angular segments. Equation ( lA.20t is solved nu- 
merically on each of these rays in the observer's frame on 151 
velocity grid points to sample v', where the transport coeffi- 
cients are pre-calculated on the grid points and later interpolated 
along the rays, whereas the profile function is always calculated 
from scratch. The numerical scheme features a variable step size 
which is controlled by comparing the results after two consecu- 
tive steps with the results obtained after one step of double size. 

A.8. Chemical heating 

By definition, exothermic chemical reactions convert chemical 
potential energies into heat, whereas endothermic reactions con- 
sume internal kinetic energy and actually cool the gas. We cal- 
culate this chemical heating/cooling rate [erg/cm 3 /s] as 



chem 



£/?(r)yf em Aff r 



(A.29) 



where r is an reaction index, R{r) is the reaction rate [l/cm 3 /s] 
and AH r [erg] is the reaction enthalpy, which is positive for 
exothermic reactions and negative for endothermic reactions 



(A.30) 



AH r = 2 Aff°(pr) - J] Aff°(ed) . 

pr ed 

AH j [erg] is the heat of formation of the chemical species in- 
volved in the reactions (pr means products, ed means educts). By 
simplification, we neglect the temperature-dependence of AHf, 
and t ake the values for the heat of formation for all species 
from dMillaretal.L[T997l) . who list AH° f [kJ/mol] at OK in their 
Table 2. 

The details of exothermic reactions are quite difficult and of- 
ten not precisely known. The excess chemical binding energy is 
seldomly released in form of kinetic energy directly (although 
there are some exceptions like, for example, dielectronic recom- 
bination which is radiationless). But in the vast majority of cases, 
the reaction will create products in some kind of electronic, vi- 
brational and rotational excited states. 

The energy temporarily stored in these excitational states can 
then be either radiated away (in which case the net heating is 
small), or subsequent collisional processes can fhermalise the 
excess energy. The ratio between these two competing processes 
depends on the critical density for collisional de-excitation, 
which depends largely on type and amount of excitation, for ex- 
ample ~ 10 15 cm -3 for permitted electronic transitions down to 
~ 10 3 cm -3 for rotational transitions. Further complications arise 
in possibly large optical depth where created photons cannot 
escape directly, scatter around and drive secondary processes, 
with a higher probability to get (partly) thermalised. To avoid 
all these complications, we have simply introduced an efficiency 
^chem m gq ( |A,29I ). We also exclude certain types of reactions 
from Eq. ( lA.29t : reactions which are known to produce or ab- 
sorb photons, cosmic ray and cosmic ray-induced reactions, X- 
ray primary and secondary reactions, reactions on grain surfaces, 
and reactions which are energetically treated in larger detail else- 
where. Some of the most important reactions are found to be 



H" + H -> H 2 + e 
H 2 + O <-> OH + H 
H 2 +H->H + H + H 
H 2 + H 2 — > H 2 + H + H 



A// r = +3.72eV 
AH r = ±0.079 eV 
A// r = -4.48 eV 
AH r = -4.48 eV, 



(A.31) 
(A.32) 
(A.33) 
(A.34) 
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most stable molecules: 
H 2 ,H 2 0, CO, CH 4 ,NH 3 



Fig. A.l. Sketch of energetics in astrochemical reaction cycles. 

where we put yj: hem = 1 for the last two "collider" reactions, as- 
suming that these highly endothermic reactions, which only oc- 
cur in hot gases, are in fact driven by the kinetic energies of the 
colliding atoms and molecules. 

In kinetic chemical equilibrium, as assumed in this paper, 
there is no net formation/destruction of any chemical species, 
but the chemistry is organised in complicated reaction cycles as 
sketched in Fig. IA.ll In the absence of cosmic rays, UV photons 
and X-ray irradiation, only the energetically most favourable 
chemical configurations, like H2,H20, CO, CH4,NH3 etc., are 
abundant. However, due to the impact of these high-energy pho- 
tons and particles, much less stable molecules, atoms and ions 
are continuously created which, via some complicated chemi- 
cal paths, come cascading down in energy again, eventually re- 
forming the abundant stable molecules. Along these paths, some 
of the chemical potential energy can be lost in form of sec- 
ondary photons. The net effect of this cycle is typically a heat- 
ing, because we have - by definition - excluded the primary UV, 
X-ray and CR reactions from Eq. (1A.291 . The chemical heat- 
ing is hence yet another way to partly thermalise the energy of 
incoming high-energy photons and particles through secondary 
exothermic reactions. 

To our surprise, the chemical heating, even with low effi- 
ciency yj: hem = 0.1, results to be an important heating process 
in protoplanetary disks, in particular at the bottom of the warm 
molecular layer where many of the near-mid IR spectral lines 
are formed, and densities are of order 10 9 - 10 11 cm -3 , which, 
as far we are aware of, has not been noticed so far in other disk 
modelling papers. 

Appendix B: The Outflow Model 

We have used the HVC [OI] 6300 A and [SII] 6731 A line lumi- 
nosities to estimate the outflow mass loss rate in two independent 
ways as suggested by Harti gan et al.l ([T995, see their Eqs. (A8) 
and (A10)). Assuming the same values for electron density (only 
necessary for the OI line) and the projected velocity of the jet as 
used by in Hartigan et al., we obtain an outflow mass loss rate 
of M out fl ow = (0.9-2) x 10~ 9 M /yr, simi lar to the mass accretion 
rat e as estimate d bv dLawson et all 120041) . 

iHollenbachl (l985) showed that it is possible to relate the 
mass outflow rate to the [OI] 63.2 yum line luminosity by assum- 
ing that each atom in the outflow passes through a single shock 
wave and that the gas with T g < 5000 K cools by radiating only 
in the [OI] 63.2 //m line. This procedure is likely to overestimate 
the resulting line flux since other cooling processes can be ac- 
tive as well, and since the jet material typically passes through 
several shock wave s close to the star . Nevertheless, we use the 
formula derived bv IHollenbachl (1 1 9 85l) to obtain a rough estimate 
of the [OI] 63.2//m line luminosity as would be expected for an 
M outflow = (0.9-2)x 10~ 9 M G /yr outflow. The result is a [OI] 63 fim 



10000 



5 1000 r 




10.0 
X [urn] 



100.0 



Fig. C.l. Dust opacities assumed in the best-fitting disk model, 
calculated with effective mixing and Mie theory according to the 
parameters listed in Table [5] Note that the absorption coefficient 
(red dotted) is larger than the scattering (black dashed) also at 
optical and near IR wavelengths, which is untypical for pure sil- 
icates. 



line luminosity equal to (0.6-1.5) times the observed value. This 
estimate suggests that a substantial part of the [OI] 63. 2 /mi line 
as seen by Herschel might originate in the outflow rather than in 
the disk. We note, however, that the peak of the HVC is shifted 
only by about -42 km/s with respect to the stellar velocity, and 
the formulae used above may not be appropriate to such unusu- 
ally low outflow/shock velocities. 

Another way to estimate the [OI] 63.2 fim e mission line flux 
from an outflow is to use the shock models of (Hartigan et al., 
2004), which predict the line ratio [OI] 6300 A/[OI] 63.2 jum 
for a variety of outflow shock parameters. The authors find 
6300/63.2 line ratios of about (0.8-2). For ETCha, the mea- 
sured HVC [OI]6300A line flux is (37-87) x 10- 18 W/m 2 , 
hence the predicted [OI]63.2yizm outflow line flux should be 
(30- 175) x 10" 18 W/m 2 , which is just consistent with the line 
flux as observed with Herschel, (30.5 + 3.2) x 10 -18 W/m 2 . If we 
include the LVC and work with the total [OI] 6300 A emission 
line flux, the results become inconsistent, i. e. the outflow model 
alone would already over-predict the measured flux. 



Appendix C: Details of the best-fitting disk model 



Figure IC.ll shows the dust opacities assumed in the best-fitting 
model, which results to be absorption-dominated and to scale 
roughly like A, except for the 10 /mi and 20 /mi silicate features. 

Figure IC.2I visualises the densities, UV radiation field 
strength, and resulting gas and dust temperatures of the best- 
fitting model. We obtain the typ ical 7d us t-pattern for passive 
disks (see e. g. iPinte et all [2009) with a shadow casted by the 
innermost disk regions and dust temperatures of the order of 
1450 K at the inner rim at 0.022 AU, and ~ 35 K at the outer disk 
edge located at 8.2 AU. The gas temperature structure r gas (r, z) 
follows the dust temperature structure in the midplane (where 
thermal accommodation dominates) but shows substantial devia- 
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Fig. C.2. Physical details of the best-fitting disk model. Upper row: total hydrogen nuclei density n<H) [cm -3 ], with overplotted 
contours for visual extinction Ay = 1, and gas temperature r gas = 1000 K, and strength of UV radiation field with respect to 
interstellar standard \- Second row: dust and gas temperature structures, TVst and 7^ (note the different scaling). Lower row: 
Most important heating and cooling processes. 
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r [AU] r [AU] 

Fig. C.3. Chemical details of the best-fitting disk model. Upper row: electron concentration n e i/"(H> and mean PAH-charge, with 
overplotted contours for ionisation parameter \og(x M<h>)- Second row: neutral carbon and CO-concentration, the CO-plot includes 
contours for dust temperature (blue) and gas temperature structure (white). Lower row: OH and H2O concentration. 
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tions at higher altitudes where the disk starts to become optically 
thin to UV radiation. 

In particular, there is a "hot finger" at relative heights z/r » 
0.3-0.6 in this model, stretching out from the inner rim to about 
4AU in radius. Here, gas temperatures of the order of 5000 K 
are achieved in this model due to PAH heating versus Fell and 
Mgll line cooling. Just below the hot atomic layer, there is a 
warm molecular layer with temperatures 300 - 1500K, heated 
by PAHs and chemical heating, balanced by CO ro-vibrational 
cooling. The outermost layers beyond 4AU are featured by an 
equilibrium between PAH and chemical heating, versus OI, OH, 
H2O and CO rotational line cooling. 

Figure IC31 shows some details about the chemical structure 
of the disk. There is a quite sudden transition of the charge of the 
PAHs from 3 to -1 where the ionisation parameter x/ n (H) is about 
100. The PAHs are mostly negatively charged then, except for 
the outer midplane where the gas almost completely neutralises. 
The neutral carbon, CO, OH and H2O concentrations are typical 
of a photon dominated region (PDR-structure). 



Appendix D: Variation of input physics 

D.1. Difference between dust and gas temperature 

A control model where the gas temperature is assumed to be 
equal to the dust temperature resulted in the line fluxes shown 
in Table ID. II We conclude that modelling the gas energy bal- 
ance, mostly leading to T gas > Ti ust in the line forming regions, 
is absolutely essential to understand the gas emission lines from 
protoplanetary disks. 

D.2. Influence of X-rays 

We have run a comparison disk model with X-ray heating 
and chemis t ry inc luded, as has recently been implemented by 
lAresu et ail d2010l) . with X-ray luminosity L x = 6 x 10 28 erg/s 
(Xmm observations bv lLopez-Santiago et al.Ll2010b . We assumed 
an X-ray emission temperature of Tx = 10 7 K and a minimum 
energy of X-ray photons of 0. 1 keV. This model does not re- 
sult in any observable changes in the calculated line fluxes. 
The modification by X-rays are only -0.5%, 1.1%, 1.4% and 
0.05% for CO / = 3 -> 2, [OI]63jum, [OI]6300A (LVC) and 
o-H2 2.122//m, respectively. Since the X-rays are attenuated by 
gas in the model, but the UV photons by dust, and our best- 
fitting model is extremely gas-rich (assumed gas/dust mass ratio 
~ 23000, see Table [5]), the X-rays do not penetrate deep enough 
to change the energy balance in the line emitting regions. 



Table D.l. Computed line fluxes [10~ 18 W/m 2 ] from a model 
where T g . ds = 7d ust is assumed in comparison to the best-fitting 
disk model. 



line 


A [pan] 


-^gas — * dust 


best model 


[OI] 'Pi -» J P 2 


63.18 


5.2 


34.5 


[OI] 'D 2 -> 3 P 2 (LVC) 


0.6300 


2.3 


69.6 


COJ = 3^2 


866.96 


0.0053 


0.014 


o-H 2 v=l->0S(l) 


2.122 


0.0073 


2.4 



D.3. Influence of viscous heating 

We have run a comparison dis k model with visc ous (gas) heating 
included, via the formula o f lFranketal.1 (119921) 

9 2 

Tvis = jPVkinO^p , (D.l) 

where p is the gas mass density, v^ n - acj H p the viscosity, a 
the viscosity parameter, cj = ^p/p the isothermal sound speed, 

Hp = cr/fikep the pressure scale height, and Qk ep = ^GM+fr 3 
the Keplerian angular velocity. Putting the viscosity parameter 
to a — 0.01, we find no significant changes in the far-IR and 
(sub-)mm lines, but modest changes in the calculated line fluxes 
in the optical and near-IR. The enhancement by viscous heating 
is 1.0%, 2.4%, 25% and 16% for CO 7 = 3 -> 2, [OI]63yum, 
[OI]6300A (LVC) and o-H 2 2.122 //m, respectively. Since the 
viscous heating scales as F v ; s oc p, but most cooling processes 
scale as A oc p 2 , the effect of the viscous heating is actually 
strongest in the low density uppermost disk regions, which is 
counter-intuitive. It renders the gas temperature in these layers 
unbound (artificially limited by 20000 K), as no implemented 
cooling process is able to balance the viscous heating according 
to Eq. ( ID. 1ft in a thin gas. We therefore refrain from discussing 
the effects of viscous heating any further in this paper. 

DA. Self-consistent disk structure 

Figure IC.4I visualises the density structure as resultant from a 
self-consistent disk model where the vertical disk stratification 
is a result of radiative transfer, chemistry, and gas energy bal- 
ance, assuming vertical hydrostatic equilibrium. We find a good 
match concerning the flaring angle (the slope in the lower left 
plot), but a generally flatter disk structure, more condensed to- 
ward the midplane. The z-dependent scale heights from the self- 
consistently calculated disk model are calculated as 



pK 2log[p(r,zVp{r,0)] ' 

where p denotes the gas pressure. In the lower left plot of 
Fig. IC.4I we have plotted two scale heights from the self- 
consistent disk model, one measured close to the midplane (at 
relative height z/r = 0.05, red dotted line) and one measured 
high above the midplane (at z/r = 0.5, blue dotted line). The dif- 
ference between these two scale heights is a natural consequence 
of our calculated gas temperature structure, with cold conditions 
in the midplane and a warm/hot disk surface. 

We observe a fair match of the scale heights at z/r = 0.5 with 
the prescribed scale-heights from our best-fitting model, but in 
the midplane the scale heights of the best-fitting model are about 
a factor of 2-3 too large. Figure IC41 also demonstrate that, con- 
sequently, we loose our SED-fit when using the self-consistent 
model. The self-consistent model intercepts less star light, re- 
sulting in a significant cooling and flux deficit in the near and 
mid IR as compared to the observations. 

D.5. Influence of non-radiative dust heating 

In our best-fitting model, we have ignored non-radiative heat- 
ing/cooling of the dust grains when determining the dust tem- 
perature structure Ti ust (r, z). However, since the gas is typically 
warmer than the dust, inelastic gas-grain collisions (thermal ac- 
commodation) lead to a collisional, i.e. non-radiative, heating 
of the dust. If we include this effect (see Eqs. (14) and (108) 
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Fig. C.4. Density structure of inner rim and impact on SED. The upper left plot shows the prescribed density structure of our 
best-fitting model. Dashed contour lines in the upper panel refer to the calculated dust temperatures. The upper right plot repre- 
sents a model with identical parameters where, however, the vertical disk structure is calculated consistently with the resultant gas 
temperatures and mean molecular weights, which results in a flatter midplane close to the star and too little near-mid IR excess 
(lower right plot). Note the "soft inner edge" (see lWoitke et al.Ll2009t) . The lower left plot compares the prescribed scale height of 
the best-fitting model (full line) with the scale heights resulting from the model with calculated vertical disk stratification. The red 
dotted line shows these results at small relative height z/> = 0.05, and the blue dotted line at z/> = 0.5. 



in lWoitke et all 120091) . we do not observe much of an effect on 
the calculated line fluxes, but the dust temperatures result to be 
slightly higher, with noticeable effects on the SED, see Fig. ID. II 

According to this model with the dust in non-radiative equi- 
librium, the temperature contrast between gas and dust is mainly 
driven by exothermic chemical reactions which are active even 
in quite deep and dense layers (see Fig. IC.2b . causing an overall 
heat transfer from gas to dust in the disk of ~4.5 x 10~ 3 L Q , i. e. 
about 5% of the stellar luminosity. It is this additional energy in- 
put that leaves the disk in form of additional mid IR continuum 
photons, causing the depicted variations in Fig. ID. II 

The effect of non-radiative dust heating on the SED is similar 
to increasing the scale heights. We have performed an additional 
run of the evolutionary strategy with enabled non-radiative dust 
heating. This run did not entirely converge. The final parame- 



ter set had a gas mass of 1.2 x 1CT 3 M , a scale height of only 
0.007 AU, i.e. a reduction of ~ 35% with respect to our best- 
fitting model. Unfortunately this is a slow and unstable option, 
because there is an additional outer iteration necessary between 
gas and dust temperature determination, to achieve consistent re- 
sults, which requires about 3 times more computational time. 

D.6. Influence of treatment of H 2 -formation 

The formation of H2 on dust grain surfaces is one of the most 
important first steps to initiate a rich molecular chemistry. It 
has profound effects also on other abundances, for instance the 
C + /C/CO transition (because of the mutual H2/C shielding) and 
on the formation of OH and H2O. Yet its rate is still rather un- 
certain. Our default choice is to calculate this key chemical pro- 
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Fig. D.l. Comparison between predicted SEDs of two models 
with and without non-radiative dust heating. The black model 
shows once more the SED of our best fitting model, without non- 
radiative dust heating. In the green model, non-radiative dust 
heating through thermal accommodation is included. All model 
parameters are identical otherwise. 



Table D.2. Calculated o-H 2 and o-H 2 line fluxes [1(T 18 W/m 2 ] 
and FWHM [km/ s] of models with different treatment of the In- 
formation on grain surfaces. 





o-H 2 2.122yum 
flux FWHM 


o-H 2 179.53 fim 
flux FWHM 


o-H 2 78.74 //m 
flux FWHM 


model A 
model B 
model C 
observed 


I. 81 30 
2.35 37 

II. 3 41 
2.5±0.1 18± 1.2 


1.41 8.2 
1.40 8.2 
1.51 8.3 
< 5.0 


11.2 8.5 

11.1 8.5 

12.2 8.6 

< 30 


model A = Sternberg & Dalgarno) d 19951). model B (best-fitting model) 
= ICazaux & Tielensl d2004h. model C =fcazaux & Tielens! (20101. The 



other calculated gas emission lines are less affected. 



cess according to (ICazaux & Tielensl 120041 called "model B" 
in Table ID. 21 ). We have run two comparison models with two 
other formulation s, one with a typically smal ler H2-formation 
rate according to (ISternberg & Dalgarnolll995[). And one w ith a 
typically larger rate according to (ICazaux & Tiele nsl, 120101). 

Co ncerning the formulation of ISternberg & Dalgarnol 
(fl99l . which is valid for standard ISM size distribution and 
dust/gas ratio only, we add a scaling factor to account for 
deviations of the total dust surface per hydrogen nucleus in 
disks as 

pform _ o v m-18», It — ( fl ") »dust/»<H> 

R m -3X10 » <H> Vr gas 1A -22„„2 • ( D ' 3 ) 



5.899 x 10- 22 cm 2 

The H2-formation rate coefficient /?^™ [1/s] needs to be mul- 
tiplied by the neutral hydrogen atom density «h to get the 
H2-formation rate in [cm -3 s -1 ]. The normalisation factor in 
Eq. ( ID. 3b results from <a 2 )isM [«dust/«<H>]iSM =5.899 x 10~ 22 cm 2 
under interstellar conditions, i. e. for a m ; n = 0.005 fim, a m ; n = 
0.25 yum, /7 = 3.5,p gl =3g/cm 3 andpd us t/Pgas = 0.01. 

Table ID. 21 shows a strong dependence of the predicted o- 
H2 2 .122^m line on the as s umed H2 -formation rate on grains. 
The [Sternberg & Dalgarno ( 1995 )-formalism result s in a fac- 



tor of 0.77 smaller and the Caza ux & Tielensl (|20 1 Ob -formalism 
in a factor of 4.8 larger line flux. This is a daunting exam- 
ple of h idden uncertainties in astroch emical modelling. We note 
that the Sternberg & Dalgarnol(ll995b -formalism gives a smaller 
FWHM 30km/s that brings this model closer to the observa- 
tions. All other calculated gas emission lines (including those of 
H2O) are less affected. 



